Fault-tolerant quantum computation with cluster states 
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The one-way quantum computing model introduced by Raussendorf and Briegel [Phys. Rev. Lett. 
86 (22), 5188-5191 (2001)] shows that it is possible to quantum compute using only a fixed entangled 
resource known as a cluster state, and adaptive single-qubit measurements. This model is the basis 
for several practical proposals for quantum computation, including a promising proposal for optical 
quantum computation based on cluster states [M. A. Nielsen, arXiv:quant-ph/0402005 accepted to 
appear in Phys. Rev. Lett.]. A significant open question is whether such proposals are scalable in the 
presence of physically realistic noise. In this paper we prove two threshold theorems which show that 
scalable fault-tolerant quantum computation may be achieved in implementations based on cluster 
states, provided the noise in the implementations is below some constant threshold value. Our 
first threshold theorem applies to a class of implementations in which entangling gates are applied 
deterministically, but with a small amount of noise. We expect this threshold to be applicable 
in a wide variety of physical systems. Our second threshold theorem is specifically adapted to 
proposals such as the optical cluster-state proposal, in which non- deterministic entangling gates 
are used. A critical technical component of our proofs is two powerful theorems which relate the 
properties of noisy unitary operations restricted to act on a subspace of state space to extensions 
of those operations acting on the entire state space. We expect these theorems to have a variety of 
applications in other areas of quantum information science. 

PACS numbers: 03.67.-a,03.67.Pp,03.67.Lx 



I. INTRODUCTION 



Overview 



One of the most surprising recent developments in 
quantum computation is the insight that quantum mea- 
surement can be used as the fundamental dynamical op- 
eration in a quantum computer 0, Q • This insight has 
significant implications for our theoretical understanding 
of how quantum computers operate, and also for the de- 
velopment of practical proposals for quantum computing. 

Historically, the first measurement-based model for 
quantum computing was the one-way quantum computer 
developed by Raussendorf and Briegel 0. A one-way 
quantum computation is performed in two stages. In 
the first stage an entangled many-qubit state known as 
the cluster state is prepared. This is a fixed entangled 
state that does not depend on the problem instance being 
solved by the computation. Indeed, when the one-way 
quantum computation is being used to simulate a quan- 
tum circuit 1 , Q shows that the identity of the cluster 
state can be made very nearly independent of the details 
of the circuit being simulated, with the only dependence 



* Electronic 



address: 



nielsen@physics.uq.edu.au 



URL: www. qinf o . org/people/nielsen/ 

t Elec tronic address: |dawson@physics.uq.edu.au| 

URL: www. physics .uq. edu. au/people/dawson/ 
1 For a review of the quantum circuit model of quantum compu- 
tation, see 0. 



being on the depth and breadth of the circuit. In the 
second stage of a one-way quantum computation a se- 
quence of single-qubit measurements is performed on the 
cluster state. These measurements are adaptive, in the 
sense that the basis in which a qubit is measured may 
depend upon the outcome of earlier measurements. Re- 
markably, this two-stage process is sufficient to simulate 
any quantum circuit whatsoever. 

More recently, an apparently quite different 
teleportation-based Q| approach to measurement-based 
quantum computation was developed by Nielsen 0. 
This approach is based on the idea now known as gate 
teleportation, introduced by Nielsen and Chuang [{| , and 
further developed by Gottesman and Chuang [g. When 
it was first introduced, teleportation-based quantum 
computation appeared to be quite different from the one- 
way quantum computer, but subsequent work 0, B @ 
has provided a unified conceptual framework in which 
both approaches may be understood. 

Although the measurement-based models of quan- 
tum computation represent an important conceptual ad- 
vance, there is an important caveat, namely, that the 
measurement-based models assume all operations are 
carried out perfectly. Since real physical systems suffer 
from noise, to make measurement-based models scalable 
we must develop techniques for combatting noise in those 
models. 

The challenge posed by noise has been met in the 
quantum circuit model of computation with the devel- 
opment of an impressive theory of fault-tolerance, pro- 
viding a large body of techniques which can be used to 
reduce the effects of noise on quantum circuits. The cul- 
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ruination of the theory of fault-tolerance is the thresh- 
old theorem, which states that for physically reasonable 
models of noise, and provided the noise is below some 
constant threshold value, it is possible to use quantum 
circuits to efficiently simulate an arbitrarily long quan- 
tum computation with arbitrary accuracy. Although the 
threshold remains to be experimentally confirmed, there 
is now general agreement that, based on our best the- 
oretical understanding of quantum noise, the threshold 
theorem solves the problem of noise in the quantum cir- 
cuit model of computation. That is, quantum noise poses 
no problem of principle for quantum computation, only 
the (very significant) practical problem of reducing noise 
levels below the threshold value. For a survey of the 
theory of fault-tolerance and the threshold theorem see 
Chapter 10 of |3j, and references therein. 

The purpose of the present paper is to develop similar 
fault-tolerant threshold results for several measurement- 
based models of quantum computation, focusing primar- 
ily on models derived from the one-way quantum com- 
puter of Raussendorf and Briegel. We refer to this en- 
tire class of models as the cluster-state model of quantum 
computation, to emphasize the crucial role played by clus- 
ter states. Note that we use the term "cluster-state model 
of quantum computation" in a rather loose sense, using 
it to denote an entire class of models based on cluster 
states. We reserve the term one-way quantum computer 
to refer to the specific model originally suggested in 0. 

Our primary motivation in studying fault-tolerance in 
the cluster-state model is to establish that the cluster- 
state model can be used as the basis for scalable practical 
proposals for quantum computation. Although several 
proposals for experimental quantum computation with 
cluster states have been made PIEJ, such proposals can- 
not be considered scalable unless fault-tolerant methods 
of implementation are developed. In this paper we de- 
velop methods for fault-tolerant computation using clus- 
ter states, methods that are applicable in a wide variety 
of practical proposals. 

Prior work on the problem of fault-tolerant compu- 
tation with cluster states has been reported in Chap- 
ter 4 of Raussendorf 's thesis This work obtained 
a threshold for a class of noise models in which Pauli 
errors occur probabilistically in a cluster-state compu- 
tation. Our threshold result applies to a more general 
noise model that is likely to be more realistic in many 
physical systems; subject to some assumptions about lo- 
cality, we allow arbitrary non-Markovian noise to occur 
in the computation, and even allow errors to occur in 
the accompanying classical computation. Thus, our work 
should be viewed as extending and complementing the 
approach taken in We note that independent work 
extending I1T| is also being undertaken by Raussendorf 
and Briegel [l2fl. 

What is it that makes proving a fault-tolerant thresh- 
old in the cluster-state model non-trivial? The obvious 
approach to proving a threshold is to take the quantum 
circuit that we want to make fault-tolerant, convert it to 



a fault-tolerant quantum circuit using the standard pre- 
scriptions, and then simulate the resulting circuit using a 
cluster-state computation. It seems physically plausible 
that noise occurring in the cluster-state model of compu- 
tation should then be corrected by the error-correcting 
properties of the original fault-tolerant circuit. 

Two difficulties obstruct this proposal. The first dif- 
ficulty is that the qubits in the cluster tend to degrade 
before they are measured. We will see that this difficulty 
can easily be overcome by building up the cluster in parts, 
so that no part of the cluster is allowed to degrade too 
much before being measured. 

The second difficulty is more serious. While it is plau- 
sible that noise in the cluster-state computation is cor- 
rected by the error-correcting properties of the original 
fault-tolerant circuit, showing this turns out to be non- 
trivial. The key is to prove that noise in a cluster- 
state simulation of a quantum circuit can be mapped 
onto equivalent noise in the quantum circuit. Provided 
that mapping has suitable properties, a threshold theo- 
rem then follows. The greater part of this paper is spent 
in constructing such a mapping. To underscore the diffi- 
culty in proving this, we mention just one interesting sub- 
tlety: we shall see that Markovian noise in a cluster-state 
computation actually maps onto non-Markovian noise in 
the corresponding quantum circuit. This and other sub- 
tleties make the task of proving a threshold technically 
challenging. 

This discussion highlights a general point worth noting. 
Our paper provides a way of taking an arbitrary quantum 
circuit and then simulating it in a fault-tolerant fashion 
in the cluster-state model of computation. We don't, 
however, provide a direct way of making a cluster-state 
computation fault- tolerant, except insofar as a cluster- 
state computation may be regarded as a special type of 
quantum circuit computation. It would be interesting to 
investigate more direct fault-tolerance constructions ap- 
plicable to an arbitrary cluster-state computation. Of 
course, for the purposes of proving the feasibility of scal- 
able quantum computation in the cluster-state model, 
the present approach is sufficient. 



B. Optical cluster states and fault tolerance 

A topic of special interest in the present paper is the 
optical cluster-state proposal for quantum computation 
suggested by Nielsen [l(|- Optical systems offer a num- 
ber of significant experimental advantages for quantum 
computing, and this proposal thus offers a very promising 
approach to experimental quantum computation. How- 
ever, the optical cluster-state proposal also differs from 
most quantum computing proposals in that it is based 
on entangling gates that only work non-deterministically. 
This non-deterministic nature poses special difficulties 
when attempting to prove a threshold for the optical 
cluster-state proposal. Following [hJ, we now briefly re- 
view some background on this proposal that will help the 
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reader understand how it fits into the present paper. 

A priori, optics offers significant advantages for the im- 
plementation of quantum computation, such as the ease 
of performing basic manipulations, and long decoherence 
times. Unfortunately, standard linear optical elements 
alone are unsuitable for quantum computation, as they 
do not enable photons to interact. This difficulty can, in 
principle, be resolved by making use of nonlinear optical 
elements 0, \^ > a t the price of requiring large nonlin- 
earities that are at present extremely difficult to achieve. 

An alternate approach was developed by Knill, 
Laflamme and Milburn (KLM) 0], who proposed us- 
ing measurement to effect entangling interactions be- 
tween optical qubits. Using this idea, KLM developed a 
scheme for scalable quantum computation based on linear 
optical elements, together with high-efficiency photode- 
tection, feedforward of measurement results, and single- 
photon generation. KLM thus showed that scalable op- 
tical quantum computation is in principle possible using 
relatively modest resources. Experimental demonstra- 
tions [T(| [HI, UiJ |20| of several of the basic elements 
of KLM have already been achieved. 

Despite these impressive successes, the obstacles to 
fully scalable quantum computation with KLM remain 
formidable. The biggest challenge is to perform a two- 
qubit entangling gate in the near-deterministic fash- 
ion required for scalable quantum computation. KLM 
propose doing this using a combination of three ideas. 
(1) Using linear optics, single-photon sources and pho- 
todetectors, non-deterministically perform an entangling 
gate. This gate fails most of the time, destroying the 
state of the computer when it does so, and so is not 
immediately suitable for quantum computation. Several 
variants of thi s g ate have already been experimentally 
demonstrated Ha IH El 13 (2) By combining the 
basic non-deterministic gate with quantum teleportation, 
a class of non-deterministic gates which are not so de- 
structive of the state of the computer is found. (3) By 
combining the gates from (2) with ideas from quantum 
error-correction, the probability of the gate succeeding 
can be improved until the gate is near-deterministic, al- 
lowing scalable quantum computation. 

The combination of these three ideas allows scalable 
quantum computation, in principle. In practice there 
are enormous obstacles to performing even a single near- 
deterministic entangling quantum gate in this fashion. 
The proposal of [l(| eliminates much of the difficulty, 
completely removing step (3), and obviating the need for 
all but the simplest versions of step (2) 2 . This is achieved 
by combining some of the simplest elements of KLM with 
the cluster-state model of quantum computation. The re- 
sulting proposal puts near-deterministic entangling quan- 
tum gates within experimental reach, and thus offers an 
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extremely promising approach to quantum computation, 
provided suitable methods for dealing with noise can be 
developed. 



C. Content of the paper 

In this paper we prove two threshold theorems for 
cluster-state quantum computation. The first is a gen- 
eral threshold theorem applicable to a variety of possible 
implementations of cluster-state computation, assuming 
that deterministic entangling gates are available. The 
second is specifically adapted to the optical cluster-state 
proposal for quantum computation. Taken together, 
these theorems show that for a wide variety of possi- 
ble physical implementations, noise poses no problem of 
principle for fully scalable cluster-state quantum compu- 
tation. 

Before describing in detail the structure of the paper, 
it is worth noting two issues that we do not fully ad- 
dress. The first of these issues is the determination of a 
numerical value for the threshold. Although we do ob- 
tain bounds on the threshold, those bounds are obtained 
through analytic methods that are far too pessimistic. 
Our philosophy is that the problem of understanding the 
threshold is best split into two parts. In the first part, 
one attempts to rigorously prove the existence of a finite 
threshold for some large class of noise models. In the 
second part, one attempts through a combination of nu- 
merical and analytic work to obtain a realistic estimate of 
the threshold for some specific and physically-motivated 
noise model. In this second part it is much more rea- 
sonable to rely on numerical evidence and heuristic rea- 
soning, since results for specific noise models can always 
be checked by computer simulation (and ultimately by 
experiment). Examples of this kind of work for quantum 
circuits may be found in [H HI |H El HI , and refer- 
ences therein. Our focus in the present paper has been 
on the first part of this program, obtaining a rigorous 
proof that a finite threshold exists. Detailed numerical 
simulation and optimization of the threshold value for 
realistic noise models is underway, and will be reported 
elsewhere p7| . 

The second issue not fully addressed in this paper 
relates to the noise model used in our analysis of the 
optical cluster-state proposal for quantum computation. 
Physically, one of the most significant sources of noise 
in any optical implementation is likely to be photon 
loss. This causes the state of the optical qubit to "leak" 
from the degrees of freedom associated with the qubit 
out into some other dimensions of the physical state 
space. In the context of threshold theorems, such noise 
is known as a "leakage error", and there are standard 
techniques for dealing with such errors in the theory of 
fault-tolerance. However, our threshold analysis for the 
cluster-state model is based on the recent threshold theo- 
rem proved by Terhal and Burkard |2S| , and that thresh- 
old does not explicitly deal with leakage errors. While it 
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seems extremely likely to us that the result of [28| can 
be patched so that leakage errors are accounted for, we 
have not worked through the analysis in detail. Rather 
than do so, in this paper we restrict ourselves to a brief 
discussion of leakage, deferring full investigation of this 
issue to a future publication. The other alternative, of 
course, would be to base our analysis on an alternative 
version of the threshold theorem. However, we will see 
that the theorem of Terhal and Burkard has considerable 
advantages for the analysis of fault-tolerance with clus- 
ter states, as the only published version of the threshold 
formulated specifically to deal with non-Markovian noise. 

The structure of the paper is as follows. We begin 
in Section [H] by defining a measure of how much noise 
occurs in a quantum information processing task. Wc 
call this measure the error strength, and prove several 
simple properties of the error strength that will be useful 
later in the paper. 

Section [H] also contains two important technical re- 
sults, which we dub the first and second unitary extension 
theorems. Roughly speaking, these results are applicable 
to situations in which two unitary operations U and V 
act in a similar fashion on a subspace S of state space. Of 
course, just because U and V act similarly on a subspace, 
it does not follow that they have similar global actions 
on state space. However, the theorems we prove guaran- 
tee that there exist unitary extensions U and V of the 
restrictions U\s and V\s, respectively, such that U and 
V have approximately equal actions everywhere on state 
space. These unitary extension theorems are critical to 
our later analysis of fault-tolerance. More generally, we 
believe that these results are of substantial interest inde- 
pendent of their application to fault-tolerance, and likely 
to find application in other areas of quantum information 
science. 

Section [H] concludes with a review of the content of 
the threshold theorem for quantum circuits. We focus 
our attention on the threshold theorem proved recently 
by Terhal and Burkard 1281, e xtending earlier work of 
Aharonov and Ben-Or [2Si l30| . The threshold of 
is unique in that it is specifically designed with non- 
Markovian noise in mind. While several of the other 
known variants of the threshold theorem can cope with 
some level of non-Markovian noise, those other variants 
are designed primarily with the case of Markovian noise 
in mind. This is important for us as we will see that 
non-Markovian noise arises naturally in the analysis of 
fault-tolerant cluster-state quantum computation, even 
if the actual physical noise occurring in the cluster-state 
computation is Markovian. 

In Section IIHI we describe the cluster-state model of 
quantum computation. Rather than providing a detailed 
proof of how the model works (which is available else- 
where) we describe the model through some simple ex- 
amples. We also discuss ways of alleviating one of the 
key difficulties that arises when attempting to perform 
fault-tolerant quantum computation with cluster states, 
the tendency of qubits in the cluster to degrade before 



they are measured. We conclude with a brief review of 
how the cluster-state model of quantum computation can 
be combined with the ideas of KLM to obtain a scheme 
for optical quantum computation. 

Section HVI is the heart of the paper, stating and prov- 
ing our first threshold theorem for cluster-state compu- 
tation. More precisely, what we prove is that it is pos- 
sible to simulate an arbitrary quantum circuit using a 
noisy cluster-state computation, with arbitrary accuracy 
and only a small overhead in the physical resources re- 
quired, subject to reasonable constraints on the physical 
resources available, and on the noise afflicting the imple- 
mentation. 

The sequence of ideas used to prove this threshold the- 
orem is conceptually rather simple. Step one is to trans- 
late the quantum circuit that we want to simulate into 
a fault-tolerant circuit, using the standard prescriptions 
for making a circuit fault-tolerant. At this stage we as- 
sume there is no noise in the computation. Step two is 
to translate the fault-tolerant circuit into a cluster-state 
computation, again using standard prescriptions. Step 
three is to carefully specify a procedure for physically 
implementing the cluster-state computation, a procedure 
that avoids the degradation of parts of the cluster that 
was mentioned above. The idea is to build up only part 
of the cluster at any given time, adding extra qubits into 
the cluster as required. Up until this step we assume 
that all operations are perfect. Step four is to introduce 
noise into the description of the cluster-state computa- 
tion, as would occur in an actual implementation. The 
most complex and critical step, step five, is to show that 
the noisy cluster-state computation is equivalent to the 
original fault-tolerant circuit, with some noise added into 
the circuit. That is, we want to map noise in the cluster- 
state computation back onto equivalent noise in the orig- 
inal fault-tolerant circuit. We will see that this mapping 
has the property that provided the noise in the cluster- 
state model is of an appropriate form, and not too strong, 
it is equivalent to noise in the original fault-tolerant cir- 
cuit which is only slightly stronger, and which is of a form 
which can be suppressed by the usual fault-tolerance con- 
structions for quantum circuits. This mapping of noise 
models thus enables us to infer a threshold theorem for 
noisy cluster-state computation. 

Section ]V\ extends these ideas to the optical cluster- 
state proposal for quantum computation. The reason 
the results of Section llVl cannot be immediately applied is 
that the entangling gates used in the optical cluster-state 
proposal are non-deterministic. We resolve this problem 
by devising an approach to computation in which the 
non-deterministic entangling gates can be treated as de- 
terministic entangling gates, subject to a small amount 
of additional noise. This enables us to map noise in the 
optical cluster-state proposal into equivalent noise in the 
deterministic cluster-state model, and then use the re- 
sult of Section llVl to infer a threshold theorem for optical 
cluster-state quantum computation. 

Section IVII concludes the paper with a summary of 
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results, and a discussion of the outlook for further devel- 
opments. 

II. NOISE AND FAULT-TOLERANT 
QUANTUM CIRCUITS 

In order to prove a threshold for cluster-state computa- 
tion we first need a way of describing quantum noise, and 
quantifying its effects. In Subsection III Al we introduce 
a measure that quantifies the effects of quantum noise, 
describe some properties of that measure, and describe 
the unitary extension theorems. In Subsection lll Bl we re- 
view the threshold theorem for quantum circuits proved 
by Terhal and Burkard [28| . 

A. Error strength 

Suppose we have a quantum system, Q, and we wish to 
implement a unitary operation Uq on that system. Un- 
fortunately, the system is not completely isolated from 
its environment, E, and thus the true evolution of the 
system will be described by some unitary evolution Vqe 
acting on both the system and the environment. We de- 
fine the error strength or noise strength of this operation 
by: 

A Q .. E (U Ql V QE ) = mm\\V QE -U Q ®U E \\, (1) 

Ue 

where the minimum is over all unitary operations Ue on 
the system E, and the norm is the usual operator norm. 

We will use the error strength Aq : e as our principal 
measure of noise in the implementation of quantum com- 
putation, whether it be by cluster-state methods, or by 
quantum circuits. Although Aq : e has been defined only 
when the ideal operation Uq is unitary, we'll see later 
that it can also be used to understand noise in opera- 
tions that may not be unitary, such as measurement and 
state preparation. Our main reason for using the mea- 
sure Aq : e is that it is the same measure that Terhal 
and Burkard use in their study of fault-tolerance [28| for 
non-Markovian noise models. We now describe several 
properties of Aq : e that arc useful later in the paper. 

Proposition 1 ( Chaining property) Let Uq, . . . , Uq 

be unitary operations on the system Q, and Vq E , . . . , Vq E 
be unitary operations on the combined system QE. Then 

m 

Aq.. e (U q ...U%,V q ]e---V$ 1 e) < Y, A Q^ U Q' V Qe)- 

3 = 1 

(2) 

The chaining property tells us that the total error 
strength of a sequence of imperfect operations is less than 
or equal to the sum of the individual error strengths. 
We note that this proposition and its proof is similar to 
Lemma 1 in Section 1.3 of |28(. 



Proof: Choose U° E so that A Q:E (U Q , V^ e ) = \\V^ E - 

U° Q ®U E \\, and define A j = V^ E -U 3 Q ^U E . A straightfor- 
ward induction on m can be used to establish the formula 

v% E ...v$ E = (U Q ...U Q n )®(U E ...U%) 

m 

+ Y / (U Q ---U Q - 1 ®U E ...U E - 1 ) x 

a- 1^.1 ...\ ;;>,.. (3) 

The result follows by subtracting (Uq . . . Uq) (g> 
{U E ...U^) from both sides of Eq. 0, and applying 
the triangle inequality. 
QED 

Later in the paper we will show that a noisy cluster- 
state computation is equivalent to a noisy quantum cir- 
cuit computation. A technique we'll use in the proof of 
this fact is to change the set of systems considered to be 
part of the environment. The following two propositions 
help us understand the behaviour of the error strength 
when the environment is changed in this way. 

Proposition 2 Let A, B and C be three quantum sys- 
tems, and let Ua,Ub,Vabc be unitary operations acting 
on systems indicated by the respective subscripts. Then 

A A :bc(Ua, Vabc) < Aab-.c(Ua ® U B , Vabc)- (4) 

Proof: The proof is immediate from the definition of 
Aq-.e and the fact that the set of unitary matrices Use 
on BC is a superset of the set of unitary matrices of the 
form Ub <£> Uc , where Ub is the (fixed) given matrix: 

Aa:bc{Ua,V A bc) = min 1 1 Vabc - Ua ® Ubc 1 1 (5) 

Ubc 

< min \\Vabc -U a ®U b ® U c \\ 

Uc 

(6) 

= Aab:c{U a ®U b ,Vabc)- (7) 

QED 

Proposition 3 Let A, B and C be three quantum sys- 
tems, and let Ua,VaBiVc be unitary operations acting 
on systems indicated by the respective subscripts. Then 

A A :bc(Ua,Vab ®V c ) < A A :b(U a ,Vab). (8) 

Proof: Similarly to the proof of the previous proposi- 
tion, the proof is immediate from the definitions and the 
fact that the set of unitary matrices Ubc on BC is a su- 
perset of the set of unitary matrices of the form Ub <%> Vc, 
where Vc is the (fixed) given matrix: 

A A :Bc(Ua,Vab®Vc) 
= min \\Vab ®V c -U a ® U B c\\ (9) 

Ubc 

< min ||Vab ® V C -U A ®U B ®V C || (10) 

Ub 

= min || Va b -U a ®U b \ | (11) 

Ub 

= Aa:b(U Ai Vab)- (12) 
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QED 

The next proposition helps in commuting noisy oper- 
ations past one another. 

Proposition 4 Let Uq and Vq be commuting unitary 
operations on a quantum system Q. Let Uqe and Vqe 
be noisy versions of these operations involving also an 
environment E. Then there exist unitaries Uqe and Vqe 
such that (a) U qe Vqe = Vq E Uq E ; Aq. e (Uq,Uq E ) < 
Aq.. e (Vq,Vqe); and Aq :E (Vq,Vq E ) < Aq.. e (Uq,U qe ). 

This proposition tells us that if Uq and Vq commute, 
then applying a noisy version of Uq followed by a noisy 
version of Vq is equivalent to applying a noisy version of 
Vq followed by a noisy version of Uq. Furthermore, the 
noise strengths in the new versions of Uq and Vq are no 
worse than in the originals, except for interchanging the 
role of Uq and Vq . 

Proof: Using the definition of Aq :E we may choose 
unitaries U E and V E , and matrices Ajj and Ay, such 
that 

Uqe = (U Q ®U E )(I + Au); (13) 
\\Au\\ = A Q:E (Uq,U qe ); (14) 

Vqe = (I + A v )(Vq®V e ); (15) 
||Av|| = Aq.. e {Vq,V Q e). (16) 

With these choices, L + Ajj and I + Ay are easily verified 
to be unitary operations on QE. We see that 

VqeUqe = {I + A V )(V Q ®V E )(U Q ®U E )(I + Au) 

(17) 

= (I + A v )(Uq®V e )(Vq®U e )(I + Au), 

(18) 

where we used the commutativity of Uq and Vq in the 
second line. The proof is completed by choosing Uqe = 
(I + A V )(U Q ® V E ) and V QE = (V Q » U E )(I + A v ). 
QED 

To prove our threshold theorems for cluster-state com- 
putation, we need two other theorems, which we call 
the first and second unitary extension theorems. These 
results are not phrased directly in terms of the error 
strength Aq. e , but we shall see later in the paper that 
these theorems have significant implications for the error 
strength. 

The first unitary extension theorem may be motivated 
by the following problem. Suppose V is a noisy unitary 
operation approximating a noiseless unitary operation U. 
(Note that U and V here act on the same state space; 
U might correspond to Uq (£> Ue in our earlier notation, 
with V corresponding to Vqe-) For some physical reason 
we are only interested in the action of U and V on some 
subspace S of the total Hilbert space. That is, we know 
on physical grounds that all inputs to the operations are 
constrained to be in that subspace. Furthermore, there 
is another unitary operation U which acts identically to 



U on the subspace S. A natural question is whether we 
can find a unitary operation V which acts identically to 
V on the subspace S, and so that V approximates U at 
least as well as V approximates U. 

Remarkably, such an extension V always exists, and 
the proof of the first unitary extension theorem shows 
how it may be constructed. We believe this theorem has 
considerable independent interest in its own right, quite 
apart from the applications later in this paper to fault- 
tolerant computation with cluster states. 

Theorem 1 (First unitary extension theorem) Let 

U, U and V be unitaries acting on a Hilbert space T. Sup- 
pose S is a subspace of T such that U and U have the 
same action on S, i.e., U\s — U\s- (Note that we do not 
assume that U and U leave the subspace S invariant, so 
U\s and U\s should be considered as maps from S into 
T .) Then there exists a unitary extension V of V\s to 
the entire space T such that 

\\V- U\\ < \\V- U\\. (19) 

The proof of this theorem is given in Appendix 1X1 
The second unitary extension theorem answers a ques- 
tion similar in spirit, but not identical, to the question 
answered by the first unitary extension theorem. Let U 
and V be unitary operations acting on a vector space 
T, with a subspace S. Suppose U\s and V\s are close, 
i.e., \\U\s — V\s\\ is small. Can we argue that there ex- 
ists a unitary operation V extending V\s, and such that 
\\U — V\\ is also small? The second unitary extension 
theorem shows that this is always true. 

Theorem 2 (Second unitary extension theorem) 

Let U and V be unitary operation acting on a (finite- 
dimensional) inner product space T. Suppose S is a sub- 
space ofT. Then there exists a unitary operation V such 
that V\s = V\s and 

\\U-V\\<2\\U\ s -V\s\\. (20) 

The proof of the second unitary extension theorem is 
given in Appendix^ Note that this theorem may easily 
be restated in the language of isometries, if that is more 
to one's taste. It is also worth noting that the second 
unitary extension theorem implies a weaker version of 
the first unitary extension theorem. In the notation of 
the first theorem, the second theorem implies that there 
exists a unitary extension V of V|s such that \\V — U\\ < 
2\\U\ s -V\s\\<2\\U-Vl 

B. Fault-tolerance in the quantum circuit model 

The threshold for cluster-state computation proved in 
this paper is based on the threshold for quantum circuits 
proved by Terhal and Burkard 28]. In this subsection 
we review Terhal and Burkard's result. We begin with a 
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description of the assumptions they make about quantum 
circuits, including the noise model, before stating their 
main theorem. Note that the noise model used by Terhal 
and Burkard is the basis for our noise model for cluster- 
state computation, described in Subsection IIV Al 

Terhal and Burkard split the total system up into three 
types of subsystem. First, there are register qubits, which 
can be controlled and used for computation. These qubits 
are present through the entirety of the computation. Sec- 
ond, there are ancilla qubits, which may also be controlled 
and used for the computation. The difference between 
register and ancilla qubits is that the ancillas may be 
brought into the computer partway through a computa- 
tion, used as part of the subsequent computation, and 
then discarded at some later time. The third type of 
system is the environment, which is not under control. 

The computation is represented by a sequence of uni- 
tary operations. Ideally, these operations would be ap- 
plied just to the register and ancilla qubits, but inevitably 
they involve some interaction with the uncontrolled en- 
vironment. It is this interaction which causes noise in 
the computer. We will find it convenient to assume the 
interaction with the environment is unitary; by making 
the environment sufficiently inclusive the laws of quan- 
tum mechanics ensure we may always make such an as- 
sumption. 

Within this framework, our noise model may be de- 
scribed as follows. Each qubit in the computer, whether 
a register qubit or an ancilla qubit, has associated with 
it its own environment. So, for example, if we label 
the qubits Qi,Q%, . . ., then the corresponding environ- 
ments would be labeled E\ , E2 , ■ ■ ■ ■ The key assumption 
we make about noise is that non-interacting qubits have 
non-interacting environments. More precisely, suppose 
as part of the computation we want to attempt some uni- 
tary gate U on qubit Qj. This might be the identity gate, 
representing quantum memory, or it might be a more 
complex gate, like a Hadamard or Pauli gate. In reality, 
this gate will be noisy, due to interactions with the envi- 
ronment. Our assumption is that the real noisy operation 
is a unitary evolution V acting on Qj and its environ- 
ment Ej, with the other qubits and their environments 
not affected. In a similar way, if we attempt a two-qubit 
operation U between qubits Qj and Qk, we assume that 
the real noisy evolution V may involve the qubits Qj,Qk, 
and the corresponding environments Ej,E^, but not the 
other qubits or their environments. With these assump- 
tions, we say the noise in a noisy circuit is of strength at 
most r\ if each ideal gate Uj in the circuit is approximated 
by a noisy gate Vj such that AQ-E Q (Uj,Vj) < rj, where 
Q is the qubit or qubits involved in the gate, and Eq is 
the corresponding environment or environments. 

We refer to the assumption that non-interacting qubits 
have non-interacting environments as the locality as- 



sumption for noise 3 . Physically, the motivation for the 
locality assumption is that each environment is well- 
localized in space, and that environments can only in- 
teract with one another when two qubits are brought 
together to interact in a quantum gate. 

Importantly, Terhal and Burkard do not make any 
Markovian assumption. That is, each environment can 
have an arbitrarily long memory. So, for example, we 
may perform a sequence of gates in which Q\ first in- 
teracts with Ei, which then passes information onto E% 
through a subsequent gate, then onto E3 through an- 
other gate, and finally corrupts qubit Q4, say. This is in 
contrast with many other variants of the threshold theo- 
rem, where Markovian noise is assumed, i.e., qubits are 
assumed to have independent and memoryless environ- 
ments. 

In addition to the locality assumption for noise, Terhal 
and Burkard make three important additional assump- 
tions about how quantum circuits are performed: 

1. It is possible to perform quantum gates on differ- 
ent qubits in parallel. Physically, this requirement 
is due to the fact that error-correction must con- 
stantly be performed on all the qubits, even if one 
is merely attempting to maintain them in memory. 
It is possible to prove that parallelizability is a nec- 
essary condition for a threshold theorem to apply. 

2. It is possible to initialize fresh ancilla qubits in the 
state |0) just prior to their being brought into the 
computation. Physically, this requirement is due 
to the fact that the ancillas are used as an entropy 
sink to remove noise from the computation. To be 
effective in this capacity they must start in a low- 
entropy state. It is possible to prove that the re- 
quirement for fresh ancillas is a necessary condition 
for a threshold theorem to apply [3l| . 

3. Excepting ancilla preparation, all dynamical oper- 
ations applied during the computation are unitary, 
up until the final measurement at the end of the 
computation. This is not a necessary feature of a 
threshold theorem, but is a feature of the threshold 
of Terhal and Burkard. 

The third assumption, that the computation is per- 
formed using only unitary operations, is rather incon- 
venient from our point of view, since the cluster-state 
model of quantum computation inherently involves many 
measurements performed during the computation. One 
feature of our proof is that it involves the replacement 
of measurements and classical feedforward by equivalent 



Terhal and Burkard consider even more general noise models, 
which may be of interest in certain circumstances. However, the 
locality assumption is sufficiently strong to cover a very wide 
class of physically interesting noise models, and so we restrict 
attention to noise models satisfying this assumption. 
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unitary operations. The reason we use the all-unitary 
model is that we need a threshold theorem which al- 
lows non-Mar kovian noise, and at present this means 
using Terhal and Burkard's all-unitary model. Future 
improvements to the threshold theorem for cluster states 
may come by developing threshold results for quantum 
circuits which allow both non-Markovian noise and mea- 
surement during the computation. 

To conclude preparation for the statement of the 
threshold theorem we need a few final items of notation 
and nomenclature. It will be convenient to assume that 
each unitary operation performed during the computa- 
tion takes the same amount of time, At, and so the cir- 
cuit may be written as a sequence of unitary operations 
performed at times t = 0, t = At, t = 2 At, and so on. We 
define a location in the circuit to be specified by a triple 
(kAt, U, Q) consisting of the time kAt at which the gate 
U is performed on a qubit or ordered pair of qubits, Q. 
Our measure of the total size of the circuit is the total 
number of locations in that circuit. Note that it is im- 
portant to count locations at which the identity gate is 
applied to a qubit. 

Computation is concluded by measuring the computer 
in the computational basis to produce a probability dis- 
tribution p. The goal of fault-tolerance is to take a perfect 
quantum circuit which outputs a probability distribution 
p and to construct a fault-tolerant quantum circuit that 
may be subject to noise, but nonetheless outputs a proba- 
bility distribution p' which is close to p in some suitably 
defined sense. As the measure of closeness we use the 
Kolmogorov distance, \\p — p'\\i = 1^ \p{ x ) ~p'( x )\- 

Theo rem 3 (Threshold theorem for quantum cir- 
cuits [2b]) There exists a constant threshold 77th > for 
quantum circuit computation with the following property. 
Let n be the number of locations in a perfect quantum cir- 
cuit, C, which outputs the probability distribution p. Let 
e > 0. We can efficiently construct a noisy quantum cir- 
cuit, C , with a total number of locations npolylog(n 2 / e) , 
and such that if the noise in C satisfies the locality as- 
sumption and is of strength at most r/ t h then the output 
distribution p' from C satisfies \\p — p'\\i < e. 

Terhal and Burkard's construction of C is based on a 
particular type of quantum error-correcting code dubbed 
a computation code by Aharonov and Ben-Or (definition 
number 15 in 30] ). As a consequence of this construc- 
tion, the circuit C is built up out of a special restricted 
class of quantum gates, gates that can be implemented 
in a fault-tolerant manner. For example, it is possible 
to construct C using just operations from the follow- 
ing restricted set: preparation of qubits in the state |0); 
the identity gate, i.e., quantum memory; H (Hadamard) 
gates, Z n /4 and Z^/g gates, where Zg is the rotation of 
a single qubit by 9 about the z axis of the Bloch sphere; 
and controlled-NOT gates. As noted earlier, C does not 
include any measurement or classical processing of data, 
except at the output; all dynamical operations are fully 
unitary. 



For our purposes in this paper it is convenient to re- 
place C with an equivalent circuit, C" , built up from a 
different set of basic operations. We make this replace- 
ment in two stages. The first stage is to replace the oper- 
ations in the circuit C by operations from the following 
set: preparations of a qubit in the state |+); the identity 
gate; gates of the form X a Zp; and the controlled- Z gate, 
which we shall call CPHASE. That this can always be 
done follows from well-known quantum circuit identities. 
We call the resulting circuit C" . 

The second stage is to replace the operations in C" by 
operations from the following set: preparations of a qubit 
in the state |+); gates of the form HZ a ; and the gate 
(H (g) H) CPHASE. We refer to this as the canonical set of 
allowed operations. To see that this can be done requires 
a little care, due to the absence of the identity gate from 
the canonical set. The trick is to simulate each gate in C" 
by two gates from the canonical set, as follows: I — ► HH; 
XaZp — > HZaHZp; CPHASE — > (H (g) H) (if ® H) CPHASE. 

We call the circuit that results when these substitu- 
tions are made C" , and refer to it as the canonical form 
of the fault-tolerant circuit C. It is clear on physical 
grounds that the canonical form also satisfies the thresh- 
old theorem. Alternately, a rigorous proof of this fact 
follows from the chaining property for error strength, 
Proposition ^ The essential idea of the proof may be 
illustrated by example: suppose an identity gate in the 
original circuit C has been replaced by two consecutive 
H gates in the canonical circuit C". Provided the H 
gates both suffer from noise of strength less than r?th/2, 
Proposition H ensures that their product is equivalent to 
doing the identity gate with error strength at most 77th- 
Thus, while the threshold n"^ for C" may be somewhat 
reduced from the threshold for C", 77th) it is reduced at 
most by some constant factor. 

Summing up, we have the following restatement of the 
threshold theorem in the form that will be used for our 
analysis of fault-tolerant cluster-state quantum compu- 
tation. 



Theorem 4 (Threshold theorem for quantum 
computation, with circuits in canonical form) 

There exists a constant threshold r/th > for quantum 
circuit computation with the following property. Let n be 
the number of locations in a perfect quantum circuit, C , 
which outputs the probability distribution p. Let e > 0. 
We can efficiently construct a noisy quantum circuit, 
C , using only operations from the canonical set of 
operations (preparation of a qubit in the state |+); gates 
of the form HZ a ; and the gate (H ® H ) CPHASE ), with 
a total number of locations npolylog(n 2 / e) , and such that 
if the noise in C satisfies the locality assumption and 
is of strength at most r)th then the output distribution p' 
from C satisfies \\p — p'\\ \ < e. 
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III. CLUSTER-STATE QUANTUM 
COMPUTATION 

In this section we describe how cluster-state quantum 
computation works. Subsec. IlII Al gives a basic descrip- 
tion of the model, and introduces language useful in the 
later analysis of fault-tolerance. Subsec. fill Bl describes 
how cluster-state computation can be realized in optics. 
All proofs are omitted, and the reader is referred instead 
to Q, or to the leisurely pedagogical account in [32| . 
Note that our account barely scratches the surface of the 
work that has been done on cluster-state computation: 
the interested reader should also consult 0> H, B EE E3 
for further work on the cluster-state model of quan- 
tum computation; further work on other measurement- 
based models of quantum computation may be found 

in p, m El H Hi ES E3 



A. Introduction to cluster-state quantum 
computation 




FIG. 1: A simple cluster state. Note that each circle repre- 
sents a single qubit in the cluster. 

The basic element of the cluster-state model is the clus- 
ter state, an entangled network of qubits 4 . An example 
of a cluster state is represented in Figure ^ Each circle 
represents a single qubit in the cluster. We may define 
the cluster state as being the result of the following two- 
part preparation procedure: first, prepare each qubit in 
the state |+) = (|0) + |l))/\/2, and then apply CPHASE 
gates between any two qubits joined by a line. Since 
the CPHASE gates commute with one another, it does not 
matter in what order they are applied. Note that this 
is merely a convenient way of defining the cluster state, 
and there is no requirement that it be prepared in this 
way. 

Given the cluster state, a cluster-state computation is 
simply a procedure for measuring some subset of qubits 
in the cluster, using single-qubit measurements and feed- 
forward of the measurement results to control the bases 
in which later qubits are measured. The output of the 



4 The states we call cluster states are in fact a generalization of 
the cluster state used in [lj. These generalized states have been 
called graph states elsewhere; we prefer to use the more elegant 
term cluster state to refer to all the states in this class. 



computation is the joint state of whatever qubits remain 
unmeasured at the end of the computation. 
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FIG. 2: A two-qubit quantum circuit. Without loss of gener- 
ality we may assume the computation starts with each qubit 
in the |+) = (|0) + |l))/\/2 state, since single-qubit gates 
can be prepended to the circuit if we wish to start in some 
other state. The two-qubit gate is a CPHASE gate. The boxes 
are single-qubit gates of the form HZ a , where H and Z a are 
as defined in Subsection III Bl Note that by composing three 
of these gates we can obtain an arbitrary single-qubit gate. 
Thus, gates of the form HZ a , together with CPHASE gates, 
are universal for quantum computation. 

Remarkably, this procedure can be used to simulate an 
arbitrary quantum circuit. Indeed, the cluster state of 
Figure n was specifically chosen in order to simulate the 
circuit in Figure |21 In Figure [3] we illustrate visually how 
the cluster state of Figure n can be used to simulate the 
circuit in Figure [2J Each qubit in the quantum circuit 
is replaced by a horizontal line of qubits in the cluster 
state. Different horizontal qubits in the cluster represent 
the original qubit at different times, with the progress 
of time being from left to right. Each single-qubit gate 
HZ a in the quantum circuit is replaced by a single qubit 
in the cluster state. CPHASE gates in the original circuit 
are simulated using a vertical "bridge" connecting the 
appropriate qubits. The cluster-state computation itself 
is carried out by performing a series of measurements in 
the time order and measurement bases indicated in the 
caption to Figure [3J The final output of the cluster-state 
computation \ip) is related to the output of the quan- 
tum circuit |0) by \ip) — cr|0), where a is a product of 
Pauli matrices that is an easy-to-compute function of the 
measurement outcomes obtained during the cluster-state 
computation. 

Although the example we have described involves a 
specific quantum circuit, general quantum circuits can 
be given a cluster-state simulation along similar lines. 
Note that we have not explicitly explained how the mea- 
surement feedforward procedure works, nor the precise 
function of measurement outcomes that determines the 
Pauli correction a at the end of the computation. These 
are explained in detail in 0, E2 ; we also give an explicit 
description of the procedures used in Section llVI 

One feature of our example quantum circuit, Figure [21 
that deserves attention is the fact that it doesn't involve 
any ancilla qubits. An important feature of many quan- 
tum circuits is that they involve the preparation and dis- 
carding of ancillas during the computation. This is es- 
pecially true of circuits for quantum error-correction and 
fault-tolerant quantum computation, where the ancillas 
are used as a heat sink to remove excess entropy from 
the computer. Such ancilla preparations and removal are 
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Simulates CPHASE 
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FIG. 3: The cluster state of Figure marked to indicate how 
it is used to simulate the different elements in the circuit of 
Figure [5] Note that the labeled qubits all have labels of the 
form im U" , where n is a positive integer, and U is a unitary 
operation. The label n indicates the time order; qubits with 
the same label can be measured in either order, or simultane- 
ously. U indicates that the qubit is measured by performing 
the unitary U and then measuring in the computational ba- 
sis. Equivalently, a single-qubit measurement in the basis 
{U' |0), I/^l)} is performed. Note that except for the first 
measurements, all measurements have U of the form HZ± a ; 
the ± indicates that the value of the sign depends upon the 
outcomes of earlier measurements. The output from the com- 
putation is at the unlabeled qubits, which are not measured. 



easily simulated in the cluster-state model. Figure^lillus- 
trates a simple quantum circuit involving an ancilla that 
is prepared and later discarded. Figure El illustrates how 
this preparation and discarding can be simulated within 
the cluster-state model. 
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FIG. 4: A simple quantum circuit in which an ancilla is pre- 
pared and later discarded. 




FIG. 5: The cluster-state computation used to simulate the 
circuit in Figure 3] 



The cluster states we have described so far have all 
been embedded in two dimensions. This is for conve- 
nience only. In practice, a more complicated topology for 
the cluster may be useful in some circumstances. This 
may be achieved either by embedding the cluster in a 
higher number of dimensions, or by nonlocal connections 
between different parts of the cluster. Fault-tolerant 
quantum circuits often involve two (or more) spatial di- 
mensions, corresponding to a three-dimensional cluster- 
state computation. Note, however, that it does not follow 
that we require the use of all three spatial dimensions to 
do a cluster-state simulation of a two-dimensional quan- 
tum circuit. We will see below that it is only necessary 
to prepare a small part of the cluster at any given time, 
and this means that the cluster-state computation may 
be performed without requiring additional spatial dimen- 
sions beyond those used in the circuit being simulated. 

Despite the different possible topologies of the cluster 
state, when convenient we shall continue to discuss clus- 
ter states as though they are embedded in two dimen- 
sions. This considerably simplifies discussion (and the 
figures), and the extensions to more complicated cluster 
topologies are in all cases obvious. 

To facilitate later discussion of fault-tolerance we 
now introduce some additional nomenclature to describe 
cluster-state computations. We call a single vertical slice 
of qubits through the cluster a layer of the cluster. For 
example, in Figure |21 the two qubits with the label 4 
form a layer. A cluster-state simulation of a circuit thus 
consists of performing single-qubit measurements, layer 
by layer. We can think of a single layer in a cluster- 
state computation as similar to the instantaneous quan- 
tum state at some fixed time during a quantum circuit 
computation. 

We call a single horizontal row of qubits through the 
cluster a level. So, for example, all the qubits on the top 
line of Figure |3 represent a level. We can think of a level 
as representing the evolution of a single qubit in a quan- 
tum circuit computation. We call a cluster-state compu- 
tation involving only a single level a single-qubit cluster- 
state computation. The motivation for this nomenclature 
is that, as we have seen, such cluster-state computations 
are used to simulate single-qubit quantum circuits. Note 
that with this definition a single-qubit cluster-state com- 
putation may involve more than one physical qubit. A 
multi-qubit cluster-state computation is one that involves 
more than a single level of qubits. 

Up until now, we have described a cluster-state com- 
putation as being composed of two steps: preparation 
of the cluster state, followed by an adaptive sequence of 
single-qubit measurements. However, it is also possible 
to implement cluster-state computations in alternative 
ways. The key observation is that we can delay prepara- 
tion of some parts of the cluster until later, doing some 
of the measurements first. So, for example we could do 
a cluster-state computation via the following sequence of 
steps: 

• Prepare layers one and two of the cluster, using |+) 



11 



preparations and CPHASE gates. 

• Perform the first layer of measurements. 

• Prepare layer three, using |+) preparations and 
CPHASE gates to adjoin the third layer of qubits 
to the second layer. 

• Perform the second layer of measurements. 

• Keep alternating the steps of preparing an extra 
layer then measuring an extra layer, until the end 
of the computation. 

We call this a one-buffered implementation of cluster 
state computation, since there is always a buffer of one 
layer between the layer of qubits being measured, and 
the most recently prepared layer of qubits. We call the 
set of qubits about to be measured the current layer, and 
the layer after that the next layer. 

For fault-tolerance the one-buffered implementation of 
cluster-state computation has a great advantage over our 
original prescription in which the entire cluster is pre- 
pared first. The reason, as indicated in the introduction, 
is that if the entire cluster is prepared first, then qubits 
which are to be measured later in the computation will 
have undergone substantial degradation by the time they 
are measured, and this will unacceptably corrupt the out- 
put of the computation . 

More generally, the one-buffered implementation illus- 
trates the important point that a given cluster-state com- 
putation may have many different implementations, i.e., 
different methods for creating the cluster and performing 
the required single-qubit measurements. When proving 
fault-tolerant threshold theorems for cluster-state com- 
putation we will need to carefully specify the details of 
the implementation used. All the implementations used 
in this paper are variants on the one-buffered implemen- 
tation. 



B. Optical cluster-state quantum computing 

Our description of cluster-state computation has been 
as an abstract model of quantum computation. As de- 
scribed in the introduction the cluster-state model also 
shows great promise as the basis for experimental im- 
plementations of quantum computation in optics [Tof . 
We now briefly describe the optical implementation of 
cluster-state computation, following |10| , and some of the 
special challenges it poses for a proof of fault-tolerance. 

The proposal of ^3 is a modified version of the pro- 
posal of Knill, Laflamme and Milburn (KLM) ^| , and 
we now briefly review some of the basic elements of 



5 We thank Andrew Childs and Debbie Leung for pointing this 
fact out. 



KLM, following the review in [lj}. KLM encodes a sin- 
gle qubit in two optical modes, A and B, with logical 
qubit states \0) L = \01)ab, and \l) L = \10) ab- (Note 
that we are using the standard Bosonic occupation num- 
ber representation on the right-hand side of these defini- 
tions, not the qubit notation, so that, for example, |01)as 
indicates zero photons in mode A, and one photon in 
mode B.) State preparation is done using single-photon 
sources, while measurements in the computational ba- 
sis may be achieved using high-efficiency photodetectors. 
Such sources and detectors make heavy demands not en- 
tirely met by existing optical technology, although en- 
couraging progress on both fronts has been reported re- 
cently. Arbitrary single-qubit operations are achieved us- 
ing phase shifters and beamsplitters. 

The main difficulty in KLM is achieving near- 
deterministic entangling interactions between qubits. 
KLM use the idea of gate teleportation 0, |(| to produce 
a gate CZ n 2/( n+1 y which with probability n 2 /(n + l) 
applies a CPHASE to two input qubits, where n is any 
fixed positive integer. When the gate fails, the effect is to 
perform a measurement of those qubits in the computa- 
tional basis. Increasing values of n correspond to increas- 
ingly complicated teleportation circuits. For this rea- 
son, KLM combine these ideas with ideas from quantum 
error-correction in order to achieve a near-deterministic 
CPHASE gate, and thus complete the set required for uni- 
versal quantum computation. 

An important property of the gate CZ n 2 ^ n+1 y is that 
in the ideal case of perfect implementation, we know when 
the gate succeeds. In particular, in KLM's implementa- 
tion procedure, success of the gate is indicated by certain 
photodetectors going "click" , while failure is indicated by 
different photodetection outcomes. We call such gates 
postselected gates to indicate that whether the gate has 
succeeded is known, and can be fed forward to later parts 
of the computation. 

The advantage of the optical cluster-state proposal 
of 0| is that it only makes use of the CZi/± and CZ^/g 
gates, both of which use relatively simple configurations 
of optical elements, and avoids the use of error-correction 
in achieving a near-deterministic CPHASE gate. This re- 
sults in a greatly simplified proposal for quantum com- 
putation. 

A key observation used in the optical cluster-state pro- 
posal is an interesting general property of cluster states. 
Suppose we measure one of the cluster qubits in the com- 
putational basis, with outcome m. Then it can be shown 
that the posterior state is just a cluster state with that 
node deleted, up to a local Z m operation applied to each 
qubit neighbouring the deleted qubit. These are known 
local unitaries, whose effect may be compensated in sub- 
sequent operations, so we may effectively regard such 
a computational basis measurement as simply removing 
the qubit from the cluster. 

This is a useful observation because when the 
C^n 2 /(n+i) ! gate fails, it effects a measurement in the 
computational basis. Thus, if one attempts to add qubits 
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FIG. 6: A single-qubit circuit containing only operations from 
the canonical set. 

to a cluster using a CZ n 2/( n+1 y gate, failure of the gate 
merely results in a single qubit being removed from the 
cluster, rather than the entire cluster being destroyed. 
|10| shows that by combining this observation with a ran- 
dom walk technique, it is possible to efficiently build up 
an arbitrary cluster state using either CZ '4/9 or CZ\u 
gates. Once this is done, all the other operations in the 
cluster-state model can be done following KLM's pre- 
scription. 

IV. FAULT-TOLERANCE WITH 
DETERMINISTIC CPHASE GATES 

In this section we prove a threshold theorem for noisy 
cluster-state quantum computation. This theorem is ap- 
plicable to situations in which the cluster can be ex- 
tended during the computation using CPHASE gates that 
are noisy, but operate deterministically. In the next sec- 
tion we extend the theorem to some situations where the 
CPHASE gates operate non-deterministically, as is the case 
for optical cluster-state computation. 

Rather than attempt to invent fault-tolerant methods 
for cluster-state computation from scratch, it is natural 
to build off the existing and rather extensive body of lit- 
erature on fault-tolerant quantum circuits. As described 
in the introduction, our strategy is to consider a cluster- 
state computation that simulates a fault-tolerant quan- 
tum circuit, and then ask if the simulated fault-tolerant 
capabilities are able to correct noise in the cluster-state 
implementation. 

We therefore begin with a quantum circuit Q and, 
instead of directly translating it into the cluster-state 
model, first encode Q as a fault-tolerant circuit Tq in 
the canonical form of Theorem^] Recall that a canonical 
fault-tolerant circuit uses only preparations of qubits in 
the state |+), single-qubit gates of the form HZ a , and the 
two-qubit gate (H <g) H) CPHASE. Using the prescription 
described in Section ITTT1 it is a simple matter to translate 
Tq into a cluster-state computation, which we denote C. 

Suppose now that C is a noisy one-buffered implemen- 
tation of C. Is C equivalent to some noisy implementa- 
tion Tq' of Tq? We will show in this section that this is 
indeed the case, and moreover that the noise is of a type 
and strength that is correctable by the fault-tolerance 
built into Tq. The noisy cluster-state computation C is 
therefore a fault-tolerant simulation of the original quan- 
tum circuit Q. It is worth noting that this noise corre- 
spondence holds for any quantum circuit and its corre- 
sponding one-buffered implementation; we don't use any 
special properties of Tq in proving the noise correspon- 
dence. 



A key construction used in establishing this noise cor- 
respondence is what we call the literal quantum circuit, 
C. The literal circuit is a quantum circuit depiction of the 
operations performed during a one-buffered implementa- 
tion of a cluster-state computation. It is a literal transla- 
tion of the one-buffered implementation, and should not 
be confused with the quantum circuit Tq being simu- 
lated. As an example, consider the single-qubit quantum 
circuit depicted in Figure The corresponding cluster- 
state computation is depicted in Figure and the literal 
circuit showing the one-buffered implementation is shown 
in Figure 03 Note that although 3 qubits appear in the 
literal circuit, the quantum circuit being simulated is a 
single-qubit computation. 

The literal circuit C offers a convenient means for de- 
scribing the effects of noise in C , and for this reason 
we have gone to some trouble in Figure |S] to depict the 
correct time-ordering of events. We have, for example, 
offset the preparation of the final |+) state, since it is 
not actually prepared until later in the one-buffered im- 
plementation, and preparation at an earlier time would 
result in considerably more noise affecting the qubit. 

The noise in C is quantified by the error strength Aq-e 
of the operations appearing in the literal circuit. The 
key result of this section is that if the worst-case noise 
strength in C is 77, then the corresponding noise in the 
quantum circuit Tq satisfies the locality assumption and 
has strength at most crj, for some constant c. Provided 
ci] < 77th we conclude that the distribution p' that results 
when we measure the output of the cluster-state compu- 
tation satisfies \\p — p'\\i < e, where p is the distribution 
output from the noise-free computation. 

The section contains three parts. Subsection IIV Al 
introduces our noise model for cluster-state computa- 
tion. Subsection IIVBI proves the noise correspondence 
described above for the simplest case when the quantum 
circuit Tq and the corresponding cluster-state computa- 
tion C are single-qubit computations. All the ideas intro- 
duced in this subsection are then extended to the case of 
a multi-qubit Tq and C in Subsection IIV Cl 



A. Noise model for a one-buffered cluster-state 
computation 

The noise model appropriate to a cluster-state compu- 
tation depends critically upon the implementation pro- 
cedure used to perform the computation. The results in 




FIG. 7: A cluster-state computation implementing the single- 
qubit circuit shown in Figure [(J We will show that the noise 
in this single-qubit cluster-state computation can be mapped 
onto equivalent noise in the corresponding quantum circuit. 
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We quantify the overall strength of noise in a one- 
buffered cluster-state computation by the worst-case er- 
ror strength in any of the unitary operations, includ- 
ing the noisy quantum memory steps in preparation and 
measurement. 



FIG. 8: The literal circuit for the one-buffered implementa- 
tion of the cluster-state computation in Figure |7| 



B. Noise correspondence for single-qubit 
computations 



this section are based on the one-buffered implementa- 
tion procedure, as described in Section Ifffl Recall that 
in a one-buffered implementation the only qubits avail- 
able in the cluster at any given time are the qubits in 
the current layer, and the next layer. The computation 
is performed by repeatedly performing the following two 
steps: (1) making all the necessary measurements on the 
current layer; and (2) using CPHASE gates to add an ex- 
tra layer of qubits into the cluster. The only variation in 
this procedure comes at the very beginning of the com- 
putation, where we need to create two whole layers of 
cluster-state qubits, and at the end, where we don't need 
to add an extra layer into the cluster. 

As in the fault-tolerance results for quantum circuits, 
our results do not allow for completely arbitrary types 
of noise. Instead, we make some physically plausible as- 
sumptions about the nature of noise in the one-buffered 
implementation. The noise model we adopt allows for 
the following types of noise: 

1. Noise in unitary dynamics: We model this in a 
manner similar to the noise model for quantum 
circuits described in Subsection III Bl Each qubit 
has its own environment, and we assume that 
non-interacting qubits have non-interacting envi- 
ronments, but make no other assumptions about 
the noise. Indeed, we can use an even more general 
noise model, in which qubits at the same level in the 
cluster state are assumed to share a common envi- 
ronment, and all we assume is that non- interacting 
levels have non- interacting environments. 

2. Noise in quantum memory: Quantum memory is 
simply the (unitary) identity operation, and we 
model a noisy quantum memory step as we would 
any other noisy unitary operation. Note that this 
type of noise affects all qubits other than the cur- 
rent layer during a round of measurements. 



3. Noise in preparation of the |+) 
this is as perfect preparation of 
noisy quantum memory step. 



state: We model 
+) , followed by a 



4. Noise in measurements in the computational basis: 
We model this as a noisy quantum memory step, 
followed by a perfect measurement in the compu- 
tational basis. 



In this section we consider the simplest case of a single- 
qubit circuit Tq made up of gates HZ aj , as shown in 
Figure A cluster-state computation C simulating Tq 
is shown in Figure [7| The establishment of the noise 
correspondence is a five-step process. To help orient the 
reader, we now outline these steps. Note that the mean- 
ing of these steps may not be completely clear upon a 
first read, but hopefully will case comprehension of later 
parts of the paper. 

1. We begin with the literal circuit C depicting a one- 
buffered implementation of C. Such a circuit is 
shown in Figure |S1 

2. The literal circuit of Figure |S| does not explicitly 
contain the classical feedforward and control that 
is performed during the cluster-state computation. 
Without taking these into account, it is not possible 
to understand the effects of noise on the computa- 
tion. Thus, we expand C to explicitly include the 
classical feedforward and control. 

3. We use a series of circuit identities to transform 
the literal circuit C into an equivalent circuit that 
contains "block" operations B aj , each of which cor- 
respond directly to the action of some gate HZ a . 
in Tq. Looking ahead, the block form equivalent to 
a perfectly implemented C is depicted in Figure IT71 
with the B aj shown in Figure IT51 

4. The threshold theorem of Terhal and Burkard in- 
volves only unitary operations. Thus, the next step 
is to replace the classical elements of the blocks 
B a . with unitary quantum equivalents to obtain a 
unitary operation QB a . . Looking ahead, the corre- 
spondence with the HZ a . of Tq is explicitly shown 
in Proposition |SJ The circuit equivalent to the lit- 
eral circuit but containing unitary blocks is shown 
in Figure 1221 with the QB aj defined in Figure I2U1 

5. It can now be shown that noise in C is equivalent 
to noise within the unitary blocks QB aj . Using 
Proposition [5] and the first unitary extension theo- 
rem, Theorem ^ we show that noise of strength r\ 
in QB a . corresponds to noise of strength at most 
crj in HZ aj , for some appropriate constant c. 

Our first task in this subsection, then, is to explicitly 
insert the classical control and feedforward operations 
into the literal circuit C of Figure |H1 and to arrange that 
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FIG. 9: A circuit with the same output as that of Figure |£] 
The new two-qubit gates that have been inserted are SWAP 
gates. 



circuit into an appropriate block form. We begin by in- 
serting perfect SWAP gates into C to obtain a more com- 
pact (but equivalent) circuit which is shown in Figure [5] 
We may assume that these swap gates operate perfectly 
as they are merely a mathematical convenience. The re- 
maining operations in Figure [51 however, are real oper- 
ations and will be subject to noise. When we need to 
emphasize that a circuit contains a mix of both real and 
perfect operations we will refer to it as an imperfect cir- 
cuit, as distinct from a noisy circuit where all operations 
are subject to noise. 

There are two classical aspects of the cluster-state com- 
putation that we have not yet explicitly included in the 
literal circuit. These are the Pauli corrections intro- 
duced by the measurement, and the classical feedforward 
of measurement results to account for these corrections. 
As an example, consider the first measurement in Fig- 
ure 00 This introduces a correction X mi to the state of 
the qubit immediately below and consequently the second 
measurement is performed in the basis HZ± a2 according 
to whether mi is or 1. In general, the Pauli correction 
is given by X X Z Z where x and z are classical variables. 
Initially x and z are both zero, and they are updated 
after each measurement by the rule 



x = z + m (mod 2) 
z' = x (mod 2), 



(21) 
(22) 



where m = 0, 1 is the measurement outcome. Subsequent 
measurement is performed in the basis HZ± aj , with the 
choice determined by whether x is or 1. 

The two classical variables x and z have been intro- 
duced into the circuit in Figure ITU1 The .Error [Tpdate 
operation EU updates x and z using Equations l|21|) 
and (|22|l . and the variable x is used to control the ro- 
tations HZ± aj . An explicit definition of EU is shown in 
Figure ^2 The circuit of Figure El can be made more 
compact by introducing the notation of Figure El for the 
classical feedforward. The resulting circuit is shown in 
Figure El 

We have made several modifications to the literal cir- 
cuit £, but the output of the imperfect circuit shown in 
Figurefullis still equivalent to the output of the noisy lit- 
eral circuit in Figure |SJ To complete the construction of 
the blocks B a . referred to earlier in Step 3, we introduce 
some additional circuit identities whose effect is to com- 
pensate for the corrections X X Z Z . This will enable us to 
make exact the correspondence with the gates HZ a . in 
the quantum circuit Tq. 



|+)-t- HZ ai ^- 
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FIG. 10: A circuit with the same output as Figure|5] but with 
the classical controls explicitly drawn in. Note that x is the 
top bit, while z is the lower bit. The new circuit notation 
involving the HZ ai and HZ a2 operations indicates how the 
classical variable x is used to control the measurement basis. 
A value of x — means that HZ a is applied, while a value 
of x — 1 means that HZ~ a is applied. Note that in the case 
of HZ ai , we always have x — 0, and so HZ ai is applied, as 
expected. 
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FIG. 11: The classical error update circuit, following Equa- 
tions and 



First, at the beginning of the first block in Figure fHfl 
we prepend the gates illustrated in Figure H*4l These are 
perfect gates, and can be prepended without changing 
the output of Figure IT31 since the initial values of x and 
z are both zero. 

Second, we modify the very final block in Figure 1131 
(which is not explicitly shown) , appending operations in- 
verse to those in Figure O as illustrated in Figure El 
The reason we may do this is as follows. If all the oper- 
ations in the cluster-state computation are implemented 
perfectly, then at the end of the computation the qubit 
would be in the state X x Z z \ip), where is the output 
from the corresponding perfect quantum circuit compu- 
tation. If we then measure this state in the computa- 
tional basis, we can compensate for the error operator 
X X Z Z by appropriate post-processing of the measure- 
ment result, i.e., by adding x to the outcome of the mea- 
surement, modulo two. This process of compensating the 
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FIG. 12: A definition to make other circuits more compact. 
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FIG. 13: The output of Figure|H]is the same as the output of 
this imperfect circuit. 
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FIG. 14: The output of Figure lT3"l is unchanged if these perfect 
gates are prepended at the beginning of the computation. 

measurement results is, however, equivalent to appending 
the perfect gates illustrated in Figure El an d dropping 
the process of compensation. 

Our third and final modification is to insert the per- 
fect gates of Figurell6l between each pair of the repeating 
blocks in Figure El Another way of stating this is that 
we insert the circuit of FigureEl a t the end of every block 
in the computation, except the last block, where it has al- 
ready been inserted, and insert the circuit of Figure ll4l at 
the beginning of every block in the computation, except 
the first, where it has already been inserted. This inser- 
tion does not modify the output of the circuit, since, as is 
apparent from Figure [TBI these gates cancel one another 
out. 

With these modifications, we see that the output of the 
noisy cluster-state computation in Figure is equivalent 
to the output of the imperfect circuit illustrated in Fig- 
ure El where the operation B a is defined in Figure IT51 

This completes the construction of the repeating blocks 
B a . The circuit shown in Figure IT7I allows us a fiction 
of identifying the topmost qubit Q as a persistent data 
qubit carrying the information in a cluster-state compu- 
tation. We will see that B a effectively performs single- 
qubit gates on Q. To see how each B a corresponds to a 
gate in the quantum circuit J-q consider the circuit iden- 
tity shown in Figure HHl This identity shows that a per- 

$ T 
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FIG. 15: The effective output of Figure 1131 is unchanged if 
these perfect gates are appended at the end of the computa- 
tion. 
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FIG. 16: We insert these perfect gates between each pair of 
blocks in Figure [TBI 
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FIG. 17: The output of this imperfect circuit is the same as 
the output of the noisy cluster-state computation in Figure|H| 



feet implementation of B a is equivalent (up to a known 
global phase factor) to the effect of applying a perfect 
HZ a gate to the first qubit, initially in the state \ip). 
Intuitively, then, we would expect that the result of the 
actual imperfections in B a would be to effect an imper- 
fect HZ a gate. That is, we obtain a way of translating 
our noisy cluster-state computation into an equivalent 
noisy quantum circuit computation. 

We have introduced the identity of Figure El f° r mo- 
tivational purposes only; we omit a proof, as we prove 
a stronger result later. Interested readers may wish to 
confirm this identity by hand. What we do now is find a 
way of showing quantitatively that noise in the imperfect 
operation B a may be mapped to noise in the quantum 
circuit Tq. 

The next step in our proof, Step 4, is to replace the 
classical elements in B aj with unitary quantum equiv- 
alents to obtain fully unitary blocks, which we denote 
QB aj . The reason for doing so will become clear in the 
final step of the proof, where we use the first unitary 
extension theorem, Theorem ^ to compare an imperfect 
QB aj to a noisy HZ aj . 
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FIG. 18: The definition of the operation B a used as the basis 
for the repeating blocks in Figure IT71 
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FIG. 19: The output of the circuit on the left is identical to 
the output of the circuit on the right. For this identity we 
assume that all operations in both circuits are done perfectly 
— there is no noise. 
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FIG. 20: The definition of the operation QB a . 



The classical elements of B a are the classical vari- 
ables x, z; the error update operation EU; and the clas- 
sical controlled U± a . These have all been replaced by 
quantum equivalents to define QB a . in Figure I2UI The 
quantum error update operation QEU is defined in Fig- 
ure |^ and U± a is defined as before but with a quantum 
control. The bits carrying x and z have been replaced 
by qubits which we label X and Z . We therefore assume 
that operations performed on these qubits are noiseless, 
including the operations in QEU. 

The output of the noisy literal circuit C is thus equiv- 
alent to the output from the imperfect circuit in block 
unitary form shown in Figure 12*2*1 The qubit labeled Q is 
persistent, and we can see that the cluster-state compu- 
tation can be thought of as the successive application of 
unitaries QB aj . These unitaries require three ancilla to 
perform an operation on Q, but the following proposition 
shows that, when all operations are done perfectly, the 
effect of each QB aj on Q is identical to that of HZ aj . 

Proposition 5 The circuit identity of Figure \2!A holds, 
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FIG. 21: The definition of the quantum error update opera- 
tion, QEU. 
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FIG. 22: The output at Q in this imperfect circuit is identical 
to the output of the imperfect circuit of Figure 1171 and thus 
is equivalent to the output of the noisy cluster-state compu- 
tation of Figure |H| 



where both circuits are assumed to be perfect. All inputs 
are assumed to be arbitrary, except the fixed |+) input, 
as shown. 

The proof of the circuit identity of Figure |23 is 
straightforward, but somewhat technical. The details are 
sketched in Appendix iBl 
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FIG. 23: The circuit identity of Proposition This is a 
perfect circuit identity, i.e., all elements in both circuits are 
assumed to be performed without any noise. 

The final step in establishing the noise correspondence 
for single-qubit computation is to use Proposition and 
the first unitary extension theorem, Theorem^ to argue 
that an imperfect implementation of QB a . is equivalent 
to a noisy operation HZ aj . Following the noise model 
of Subsection llV Al let E be the environment responsible 
for the noisy operations in the imperfect implementation 
of QB a . . We denote this imperfect operation by QB' a , 
a unitary acting on QMjXZE. By assuming that the 
same environment E is reused by all the imperfect oper- 
ations QB' a . we make the most pessimistic assumption 
we can possibly make about noise in the cluster-state 
computation. In the more general multi-qubit situation 
this assumption will correspond to assuming that all the 
qubits in the same level of the cluster share the same 
environment. 

Now, suppose r\ is the maximal error strength in any 
of the noisy operations making up the imperfect QB aj . 
That is, 77 quantifies the strength of the noise in the 
cluster-state computation. Let c be the number of noisy 
operations in the imperfect QB aj . Then the chaining 
property, Proposition ^ implies that 



&QM i XZ:E(QB aj ,QB' ) < cq, 



(23) 



where QB a . is the perfect QB a - operation. Note our 
convention that in algebraic expressions we always use 
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QB' a to refer to the imperfect operation, while in the 
text we sometimes use QB aj to refer to the imperfect 
operation, provided the context is clear. It follows that 
there exists a unitary operation U E on E such that 



\QB' a .-QB aj ®U E \\ < 



cq. 



(24) 



(Note that U E may depend on j, but that dependence 
is not important in the argument that follows, and so is 
suppressed.) 

So, if QB' a . is the noisy implementation of QB a . , can 
we show that this corresponds to a noisy implementa- 
tion of HZ a . ? In Figure we see that the qubit Mj 
is always initially prepared in the state |+). This is an 
exact statement, since imperfect preparation of the state 
|+) is modeled as a perfect preparation, followed by a 
noisy quantum memory step that is absorbed into the 
imperfect operation QB a . . Define Sj to be the subspace 
of QMjXZE in which Mj is in the state |+), and the 
other systems may be in an arbitrary state. Then define 
Uj = QB aj <E>Ue and Vj = QB' a . . Define Uj as shown in 

Figure |2~H By Proposition |S] we see that Uj\ s = Uj\ s 
and so we can apply the first unitary extension theorem, 
Theorem ^ to conclude that there exists a unitary ex- 
tension Vj of Vjlsi such that 



= \\QB' a .-QB aj ®UE\ 
< cq. 

Applying Proposition [5] we see that 

^■Q:M 3 XZE{HZ aj ,Vj) < cq. 



(25) 
(26) 
(27) 



(28) 



Next, define M = Mi ® M2 ® ... to be the combina- 
tion of all the systems Mj. Applying Proposition El and 
Equation l|28|) we see that 



A Q: E>(HZ aj ,Vj)<Cr], 



(29) 



where E' = MXYE is the effective environment for the 
qubit Q, and we have extended Vj to act in the natural 
way on E' , i.e., Vj acts trivially on systems Mk such 
that k 7^ j. Because Vj\s 3 — Vj\s d we see that the noisy 
implementation Vj — QB' a of QB aj is exactly equivalent 

to a noisy implementation Vj of HZ aj of strength at most 
cq. 

To conclude, we have shown that if the one-buffered 
cluster-state computation depicted in Figure [5] is imple- 
mented with noise of strength at most 77 in each opera- 
tion, then the output of that computation is equivalent 
to the output of the noisy quantum circuit in Figure 1251 
where each operation is performed with noise of strength 
at most cq. As we have described it, c is the number 
of noisy operations in the imperfect operation QB a , i.e., 
c w 10 1 . In actual implementations it would be possi- 
ble to directly evaluate the total strength of the noise in 
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FIG. 24: The definition of the operation Uj The environment 
E is shown at the bottom, where the wire with a slash through 
it indicates an arbitrary quantum system. 



the imperfect operation QB a , resulting in a more accu- 
rate (and better, from the point of view of the threshold) 
value for c. 

A interesting feature of our argument is that even if the 
various physical operations involve only Markovian noise, 
the corresponding effective noise in the implementation 
of HZ aj is inherently non-Markovian. The reason is that 
the qubits X and Z associated with the classical variables 
x and z are part of the effective environment of Q at every 
stage of the computation, due to the necessity of feeding 
forward the measurement results. This is the reason we 
need to use the non-Markovian threshold result of Tcrhal 
and Burkard [23. 
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FIG. 25: This noisy circuit is equivalent to the noisy one- 
buffered cluster-state computation of Figure |H] 



Noise correspondence for multi-qubit 
computations 
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FIG. 26: An illustrative two-qubit quantum circuit Tq using 
canonical fault-tolerant gates. 

In this subsection we extend the ideas of the previ- 
ous subsection, explaining how noise in a multi-qubit 
one-buffered cluster-state quantum computation may be 
mapped to equivalent noise in the corresponding quan- 
tum circuit. The ideas used to do this are the same as 
were used in the proof for single-qubit cluster-state com- 
putations. As a result, we merely sketch out how the 
proof goes for a specific example, with the general proof 
following similar lines. 
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FIG. 27: A noisy two-qubit cluster-state computation. 
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FIG. 28: The output of this imperfect circuit is equivalent to 
the output of the noisy two-qubit cluster-state computation 
in Figure |53 The operations QB aj , QBp j are as denned in 
the previous subsection, while QC is denned in Figure 12^1 



The noise correspondence is established via the same 
process followed in Subsection IIV Bl We begin with a 
multi-qubit quantum circuit like the circuit Tq shown 
in Figure 1261 The cluster-state computation C simu- 
lating this circuit is shown in Figure [Z7\ As before, 
we construct a literal circuit and rearrange it into the 
block form where each block is identifiable with a gate 
in the quantum circuit J-q. We have omitted the de- 
tails of this rearrangement, and presented the final block 
form equivalent to the literal circuit in Figure 1281 The 
blocks QB aj , QBp. correspond to the single-qubit gates 
HZ ajl HZ^ j in !Fq, as in the previous subsection, and 
the effects of noise in these blocks has already been con- 
sidered. The new element is the unitary block QC, which 
is shown in detail in Figure 1291 The following proposi- 
tion shows that a perfect implementation of QC effects 
the operation {H ® i?)cPHASE on the qubits Q\ and Q 2 . 

Proposition 6 The circuit identity of Figure holds, 
where both circuits are assumed to be perfect. 

Proof: The proof of this identity uses the same tech- 
niques as the proof of Proposition|Sl and is omitted. The 
identity may easily be verified using any of the standard 
computer algebra packages. 

QED 

This proposition plays a similar role to that played 
by Proposition [S] in establishing the single-qubit noise 
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FIG. 29: The definition of the imperfect operation QC ap- 
pearing in Figure 1251 We have prepended and appended the 
same fictitious elements used in the single-qubit case in Fig- 
ures^] and ^| The gate QEU' is the two-qubit extension of 
the error update, and is shown in detail in Figure l3"0l The only 
noisy operations are the Hadamard gates, and the CPHASE 
gate between the first and the fifth qubits. 
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FIG. 30: The definition of the two-qubit quantum error up- 
date operation analogous to that of Figure |2"T1 



correspondence. In this case the noise correspondence 
will again follow from the first unitary extension theorem. 

Following the noise model of Subsection IIV Al we in- 
troduce environments E\ and E2 associated with the re- 
spective levels in the cluster-state computation of Fig- 
ure 1271 Under the locality assumption an imperfect 
QB ai , for example, is thus represented as a unitary op- 
erator QB' acting on Q\M\X\Z\E\, which results in 
an effective environment for the corresponding BIZ' of 
Ei = M\X\Z\E\. As in the single-qubit case, the noise 
strength in the operations of the form HZ a is at most 
cq, where c « 10 1 is a small constant, and 77 is the noise 
strength in the operations used to implement the cluster- 
state computation. 

An imperfect QC is similarly represented as a uni- 
tary QC' which involves both the environments E± and 
E 2 - Using the first unitary extension theorem we can 
show that this results in an effective environment E\i = 
NhX 1 Z 1 M i X 2 Z 2 E 1 E2 for the operation (If <g>.ff)cPHASE 
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FIG. 31: The circuit identity of Proposition |S| This is a 
perfect circuit identity, i.e., all elements in both circuits are 
assumed to be performed without any noise. 



in Tq. It follows that the noise strength in the (H ® 
i?)CPHASE gate is at most c'ry, where d is the total num- 
ber of noisy operations in the imperfect operation QC, 
and again is of order 10 1 . 

Summing up, suppose we perform a noisy one-buffered 
implementation of a cluster-state computation, satisfying 
the noise model described in Subsection IIV Al and with 
r\ the maximal noise strength in any operation. Then 
we can show that this noisy computation is equivalent 
to performing the corresponding quantum circuit with 
noise satisfying the locality assumption, and of strength 
at most c"rj, where c" w 10 1 is some constant. This com- 
pletes the proof that quantum circuits may be simulated 
in a fault-tolerant fashion using a one-buffered cluster- 
state computation. 

One remarkable feature of our proof is that it goes 
through unchanged even if we allow noise to occur in the 
classical computations and feedforward. It is easy to see 
that such noise simply causes an additional contribution 
to the strength of the noise in the corresponding quan- 
tum circuit, and thus a decrease in the effective threshold. 
This feature of the proof also carries over to the threshold 
theorem for optical cluster-state computation presented 
in the next section. Thus, our results show that not only 
can cluster-state computation be made resilient against 
the effects of unitary non-Markovian errors, it can even 
be made resilient against the effects of noise in the clas- 
sical parts of the computation. 



V. FAULT-TOLERANCE WITH OPTICAL 
CLUSTER STATES 

In this section we explain how the ideas of Section IIVI 
can be extended to enable fault-tolerant simulation of 
quantum circuits using optical cluster states. The main 
challenge in proving this result is the non-deterministic 
nature of the entangling gates used in optical cluster- 
state computation. We show that this challenge can be 



met by using those non-deterministic gates to add addi- 
tional pieces to the cluster in a near-deterministic fash- 
ion using what we call the dangling node implementation 
of optical cluster-state quantum computation. On those 
rare occasions when the addition to the cluster fails, one 
simply accepts failure and moves on, regarding the failure 
as a small amount of additional noise in a deterministic 
preparation of the cluster state. This allows us to apply 
the results of the Section IIVI to deduce a fault-tolerant 
threshold for optical cluster states. 



Although we phrase our results in terms of optics, there 
may be other natural contexts in which our results apply. 
The crucial fact about the optical implementation that 
we use is the existence of a postselected non-deterministic 
CPHASE gate which, when it fails, effects a measurement 
in the computational basis. The threshold results we 
apply would apply equal well to other implementations 
which share this feature. 



We begin our account in Subsection IV Al with an 
overview of the main elements of the proof. In particular, 
we explain the dangling node implementation in detail, 
and explain heuristically how noise in the dangling node 
implementation maps to noise in the quantum circuit 
model. Subsection IV HI discusses the noise model we use 
in our description of optical cluster-state computation. 
An important issue in this discussion is noise due to pho- 
ton loss. As described in the introduction to this paper, 
the noise model used by Terhal and Burkard j2£| does 
not explicitly allow for leakage errors, and thus noise due 
to photon loss is not directly addressed by our results. 
Nonetheless, we expect that with some straightforward 
modifications of the results in [2^ it will be possible to 
cope with the problem of photon loss and other forms of 
leakage error in optical cluster-state computation. Sub- 
section proves a rigorous threshold for what we call 
the two-at-a-time implementation of cluster-state com- 
putation. The two-at-a-time implementation makes use 
of deterministic CPHASE gates, and thus is not of imme- 
diate relevance for optical cluster-state computing. Its 
interest arises from the fact that it is intermediate in com- 
plexity between the one-buffered implementation studied 
in Section llVl and the full dangling node implementation 
of optical cluster-state computation. Analysing the two- 
at-a-time implementation thus provides a useful stepping 
stone on the way to understanding fault-tolerance in the 
dangling node implementation. Subsection IV Dl intro- 
duces some useful nomenclature for describing postse- 
lected quantum gates, and proves a simple lemma about 
such gates. The threshold proof for optical cluster-state 
computation is completed in Subsection IV El which ex- 
plains how noise in the dangling node implementation 
may mapped to equivalent noise in the quantum circuit 
model. 
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FIG. 32: The microcluster used in the dangling node imple- 
mentation of cluster-state computation. There are k dangling 
nodes, in general; in this example, k = 3. 



A. Overview 

The broad structure of the threshold proof for optics 
is similar to the proof for the deterministic one-buffered 
implementation discussed in Sect ion HVI We take a quan- 
tum circuit that we wish to simulate, and turn it into a 
fault-tolerant quantum circuit, using the standard pre- 
scriptions for fault-tolerance. We then convert the fault- 
tolerant circuit into an equivalent cluster-state computa- 
tion, and then specify in detail a specific implementation 
protocol for performing that computation. In this case 
that will be a dangling node implementation. We now 
describe how the dangling node implementation works. 
We will show in later subsections that noise in the dan- 
gling node implementation can be mapped onto equiva- 
lent noise in the original fault-tolerant quantum circuit. 
This enables us to deduce that the output of the noisy 
optical cluster-state computation is equivalent to the out- 
put of a noisy fault-tolerant quantum circuit computa- 
tion, and thus to obtain a threshold result. 

We begin by describing the dangling node implemen- 
tation for the simplest case of a single-qubit cluster-state 
computation. The basic idea is to to build the cluster 
up using small cluster states that we call microclusters® . 
The basic microcluster is illustrated in Figure 1321 It in- 
cludes a single base node, on the left, and k dangling 
nodes, on the right, where k is some fixed constant. 

Optically, we may prepare microclusters by sequen- 
tially performing k non-deterministic CPHASE gates. 
These may be either the non-deterministic KLM CPHASE 
gate, or one of its more efficient (but still non- 
deterministic) modern descendants. Since k is a constant, 
the expected number of operations required to form a mi- 
crocluster is also constant. In practice, there are likely 
to be much more efficient means of preparing a micro- 
cluster than this procedure. However, for the purposes 
of this paper we shall not be concerned with optimizing 
the formation of the microcluster. 

The dangling node implementation of a single-qubit 
cluster-state quantum computation works as follows. 
The first step is to prepare a single microcluster, as illus- 
trated in Figure 123 This forms the basis from which a 



A similar microcluster construction was used for different pur- 
poses in If Ol . 




FIG. 33: The cluster state which results after successfully ad- 
joining a microcluster to the initial microcluster in Figure l3~2l 
deleting the extra nodes, and applying appropriate local op- 
erations. 




FIG. 34: The cluster state which results after failing k — 1 
times to adjoin a microcluster to the initial microcluster in 
Figure 01 



larger cluster state will gradually be grown by adjoining 
extra microclusters. 

The second step of the dangling node implementation 
is to attempt to adjoin microclusters to each of the first 
k—1 dangling nodes of the first microcluster (i.e., all but 
the last of the dangling nodes), using non-deterministic 
CPHASE gates. If one of these attempts to adjoin should 
succeed, then we stop, and use computational basis mea- 
surements and single-qubit operations to ensure that we 
end up with the larger cluster state shown in Figure 1551 
To ensure that this works, it is critical that we use one of 
the special KLM non-deterministic CPHASE gates, which, 
as noted in Section ITTTl have the property that failure of 
the gate results in a computational basis measurement, 
and removal of a qubit from the cluster. If all k — 1 at- 
tempts to adjoin fail, then we simply accept the failure, 
and apply the appropriate single-qubit operation to end 
up with the larger cluster state shown in Figure [Jj] Note 
that if the non-deterministic KLM CPHASE gate fails with 
probability pf, then the probability of successfully ad- 
joining a microcluster, as in Figure E3 is while 

the probability of failing, as in Figure 1541 is Pf 1 - When 
the adjoinment succeeds, we effectively add two perfect 
extra layers of qubits to the cluster. Failure introduces 
a defect into the cluster, but we can ensure this occurs 
with small probability by choosing k to be large. 

The third step of the dangling node implementation is 
to perform the first two layers of single-qubit measure- 
ments in the cluster-state computation. We then go back 
to the second step, effectively adding another two layers 
of qubits into the cluster, and so on, through the entire 
course of the computation. The only variation comes at 
the very end of the computation, where there is no need 
to adjoin a microcluster of the form in Figure 1321 but 
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FIG. 35: The simpler two-qubit micro-cluster adjoined to the 
cluster state at the very end of a single-qubit cluster-state 
computation, in place of the microcluster in Figure 13*21 
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k dangling nodes 



FIG. 36: This microcluster is used in the dangling node im- 
plementation of multi-qubit cluster-state quantum computa- 
tions. 



instead we can add the simpler two-qubit microcluster 
illustrated in Figure 1331 

We have seen how a single-qubit cluster-state computa- 
tion may be performed in the dangling node implementa- 
tion; what about multi-qubit cluster-state computations? 
The key is to introduce a third type of microcluster, il- 
lustrated in Figure 1201 Once again, note that these mi- 
croclusters may be prepared in constant expected time 
using non-deterministic CPHASE gates. 

Consider the cluster-state computation depicted in 
Figure 1371 We may implement this computation as fol- 
lows. The first step is to prepare the first two layers of the 
cluster. This may be done using the non-deterministic 
CPHASE gate in constant time, as described earlier. 

The second step of the dangling node implementation 
is to add in two additional layers to the cluster. This 
is done in a similar fashion to the adjoinment of a mi- 
crocluster in the single-qubit case. In this case, success 
requires that we adjoin the microcluster in Figure 137)1 to 
two different dangling nodes, one on the top level of the 
computation in Figure I3T1 the other to the bottom level. 
Success in adjoinment at both levels results in the cluster 
of Figure I3"81 being prepared. A simple calculation shows 
that this occurs with probability \ — 2p k ^~ 1 +p 2k ~ 2 , which 
approaches 1 rapidly as k becomes large. If, on the other 

FIG. 37: An example cluster-state computation. We have 
omitted explicit description of the measurements performed, 
since it is the formation of the cluster that we wish to con- 
centrate on, rather than the details of the measurements. 
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FIG. 38: The cluster after successfully adding the microclus- 
ter of Figure l3l)l 



o — o 
o — o 




FIG. 39: The cluster after failing to add the microcluster of 
Figure EH 



hand, either adjoinment should fail, we simply declare 
failure overall, resulting in the cluster of Figure |2H being 
prepared. Note that if one part of the microcluster is suc- 
cessfully adjoined, but the other part fails, then it may be 
necessary to delete some nodes of the cluster using com- 
putational basis measurements and local Z operations in 
order to obtain the cluster of Figure 051 

The third step in the dangling node implementation 
of the computation in Figure 1371 is to perform the first 
two layers of single-qubit measurements. The fourth step 
is to add in the final two layers of the cluster, which is 
done using the same procedure as for single-qubit com- 
putations. The fifth and final step of the implementation 
is to do the final four layers of single-qubit measurements 
in the standard way. 

The generalization of the dangling node implementa- 
tion to arbitrary multi-qubit cluster-state computations 
follows the same lines, alternating attempts to adjoin mi- 
croclusters with two layers of single-qubit measurements. 
Each successful adjoinment of a microcluster adds two 
perfect extra layers to one or more levels of the cluster, 
while failure to adjoin correctly introduces a defect into 
the cluster, but occurs with probability 2p^ _1 — Pf k ~ 2 - 

How does noise in the dangling node implementation 
map to noise in the quantum circuit model? A heuris- 
tic argument is as follows. Roughly speaking, each gate 
in the original circuit is simulated by adjoining a micro- 
cluster and performing some single-qubit measurements. 
Adjoining a microcluster involves up to c\k 2 physical 
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operations, where c\ is some positive constant in the 
range 10°-10 2 . If each operation is performed with er- 
ror strength 77, then the total associated noise is at most 
c\k 2 rj. There is also an intrinsic failure probability, due 
to the non-determinism of the gates, which causes defects 
in the cluster. This probability is bounded by 2p k j~ 1 . Fi- 
nally, the noise contribution due to the single-qubit mea- 
surements scales as C2?y, where C2 is a constant of order 
10°. In consequence we expect that noise of strength rj 
in the dangling node implementation will map to equiv- 
alent noise of strength C\k 2 r\ + c 2 ?7 + 2p k j _1 in the quan- 
tum circuit model. It follows that if rjth is the threshold 
in the quantum circuit model then provided rj satisfies 
c\k 2 r\ + C2i] + Zp^ 1 < ?7th fault-tolerant computation 
is possible using optical cluster-state computation. The 
goal of the next four subsections is to make this heuristic 
argument rigorous. We will see that the rigorous con- 
clusions are in qualitative agreement with this heuristic 
analysis, with some minor quantitative changes. 

An issue we have glossed over in our discussion is that 
the dangling node implementation restricts the structure 
of the clusters that may be formed. In particular, con- 
nections between levels of the cluster can only be formed 
within odd numbered layers of the cluster. This has the 
effect of slightly restricting the quantum circuit opera- 
tions that may be directly simulated with such a cluster. 
Using an argument analogous to that in Subsection III Bl 
one can verify that this restricted canonical set of opera- 
tions can be used to prove a threshold theorem analogous 
to Theorem 0] We omit the details of this argument, 
which is straightforward. 



B. The noise model in optical cluster-state 
computation 

We assume that noise in optical cluster-state compu- 
tation follows essentially the same model as was intro- 
duced in Subsection IIV Al for cluster-state computation 
with deterministic CPHASE gates. In particular, we as- 
sume that noisy preparation and measurement can be 
modeled as perfect operations, accompanied by noisy 
quantum memory steps. To model noisy unitary opera- 
tions we assume that each level in the cluster has its own 
environment, and that non-interacting levels have non- 
interacting environments. The only new element that 
needs to be accounted for is when unitary operations 
are performed which involve the ancilla used during the 
non-deterministic CPHASE gate. To cope with this, we 
make the pessimistic assumption that these ancilla share 
a common environment with whichever levels of the clus- 
ter are involved in the attempted CPHASE gate. 

This noise model omits a significant possible source of 
noise, that of photon loss. The basic problem, as alluded 
to earlier, is that the noise model we are using, based 
on that in 28], does not allow leakage errors. Noise is 
assumed to arise from the interaction of a qubit with 
some environment. In reality, noise may sometimes have 



a rather different nature. In particular, a qubit is some- 
times a two-dimensional subspace of a larger physical 
state space, and noise may be due not to interaction with 
an environment, but rather to leakage of the state of the 
qubit into some other part of the state space. This is 
exactly the type of error caused by photon loss in optical 
quantum computation. 

How to deal with this type of leakage error is well- 
understood in the theory of quantum error-correction, 
and is addressed by several of the standard threshold the- 
orems. In particular, loss detection techniques are used 
in the threshold analysis accompanying the original 
KLM proposal for optical quantum computation. This 
type of error is not, however, explicitly addressed by the 
threshold theorem of |28|. Although we believe that the 
result of j2^| can likely be adapted to cope with such 
leakage errors, we have not yet performed a complete 
analysis. Such an analysis will appear in future work. In 
the meantime, our results apply to the more restricted 
model of noise without leakage. 



C. The two-at-a-time implementation of 
cluster-state quantum computation 

In Section ^3 we explained how noise in a one-buffered 
implementation of cluster-state computation may be 
mapped to noise in a quantum circuit computation. In 
this section we'll extend those results to a more complex 
implementation of cluster-state computation, which we 
call a two-at-a-time implementation. 

A two-at-a-time implementation is similar to the one- 
buffered implementation, except now qubits are added 
into the cluster two layers at a time, and the single-qubit 
measurements are performed two layers at a time. In 
particular, we assume that deterministic CPHASE gates 
are available in a two-at-a-time implementation. More 
explicitly: 

• Prepare layers one through four of the cluster, using 
|+) preparations and (deterministic) CPHASE gates. 

• Perform the first two layers of measurements. 

• Prepare layers five and six, using |+) preparations 
and CPHASE gates to adjoin the extra layers. 

• Perform the third and fourth layers of measure- 
ments. 

• Keep alternating the preparation of two extra lay- 
ers with the measurement of two extra layers, until 
the end of the computation. 

The literal circuit for a two-at-a-time implementation of 
a single-qubit cluster-state computation is shown in Fig- 
ure I4UI In this figure we have explicitly drawn in the 
noisy quantum memory steps, using little crosses to indi- 
cate when a quantum memory step is being performed. 
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FIG. 40: The literal circuit for a two-at-a-time implemen- 
tation of a single-qubit cluster-state computation, with the 
noisy quantum memory steps explicitly indicated by crosses. 
For simplicity we have combined the single-qubit rotations 
and computational basis measurements into a single-qubit 
measurement in some other basis, not explicitly specified. 
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FIG. 41: The output of this circuit is equivalent to the output 
of the literal circuit for the two-at-a-time implementation in 
Figure l40l 



Including these explicitly makes it easier to map this im- 
plementation onto an equivalent one-buffered implemen- 
tation. 

Commuting the operations in Figure I4UI forward in 
time, we see that that the output of Figure 0U] is equiv- 
alent to the output of the circuit in Figure [31] Note 
that doing this commutation requires the use of Propo- 
sition since noisy operations on different qubits that 
are on the same level of the cluster may potentially in- 
volve the same environment, and so may not commute. 
Fortunately, Proposition^] allows this commutation to be 
performed without increasing the strength of the under- 
lying noise. 

Inspection of Figure 1411 reveals that this circuit may 
be regarded as the literal circuit for a one-buffered im- 
plementation of a cluster-state computation, and thus is 
equivalent to a noisy quantum circuit, using the results 
of Section IIVI A similar argument can be used to map 
noise in a multi-qubit two-at-a-time implementation of a 
cluster-state computation to equivalent noise in a quan- 
tum circuit. We omit the details of the argument, which 
is a straightforward extension of the single-qubit case. 



D. Postselected gates 

In order to map noise in optical cluster-state compu- 
tation into the quantum circuit model, we need a simple 
lemma about postselected quantum gates. In this sub- 
section we provide a formal definition of what we call a 



unitary postselected gate, and prove the required lemma. 
The unitary postselected gates discussed here differ in an 
important respect from the postselected gates in the dis- 
cussion of optical quantum computing in Section IlIII In 
the optical gates, success of the quantum gate is condi- 
tional on some measurement outcome occurring. In the 
present scenario the postselected gate is all-unitary, i.e., 
no measurement is involved. The connection between the 
two types of postselection is made by replacing the mea- 
surement by an equivalent unitary process; we will see an 
explicit example of how this works in the next subsection. 

Let U be a unitary gate acting on two registers, labeled 
A and B. B is initially in some fixed state |/3), while A 
may be in an arbitrary state \ip). Let V be a unitary 
gate acting on register A alone. U is said to be a unitary 
postselected gate that implements V with probability p if 
for all 



urn® = y/pvmp) + y^~3m\n, m 

where \j3') is some fixed state, and \j3") is orthonormal 
to 

Lemma 1 Let U be a unitary postselected gate imple- 
menting V with probability p, on registers A and B , and 
when |/3) is input to the register B. Let S be a subspace 
of the state space for register A, and let T be the subspace 
of the total state space for AB that contains states of the 
form \ip)\/3), where \ip) € S. Then there exists a unitary 
W acting on register B such that 



\\U\T-(V®W)\ T \\ = y/2(l-y/p). (31) 

Proof: By definition, for all states \ip) in register A 
we have 



urn® = vWMjSO + V^pWW)- (32) 

Let W be a unitary operator taking \j3) to \ j3'). Then 

{u-v®wm\p) = (VF-i)F|V)|/3') 

+ ^T~pW)\p"). (33) 

Evaluating the norm, and restricting to normalized states 
I'f/')!/?) which are in T, we obtain the result. 
QED 



E. How noise in optical cluster-state computation 
maps to noise in the quantum circuit model 

In this subsection we explain how noise in the dan- 
gling node implementation of optical cluster-state com- 
putation is mapped to equivalent noise in the two-at-a- 
time implementation described in Subsection IV CI The 
results in Subsection IV CI may then be used to map that 
noise to equivalent noise in the original quantum circuit, 
which enables us to complete the threshold proof. Our 
main focus here is on noise in single-qubit cluster-state 
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FIG. 42: The literal circuit for a dangling node implemen- 
tation of a single-qubit cluster-state computation. We have 
chosen to use an implementation with k = 3 dangling nodes; 
the circuit for other values of k follows similar lines. 



of attempting to add two extra layers of qubits to the 
cluster of Figure B2"l results, with high probability, in the 
cluster of Figure OS The identity fails to be exact be- 
cause (a) the attempt to add extra layers sometimes fails, 
resulting (even in the ideal case of perfect operations) in 
the cluster of Figure 1341 and (b) the physical operations 
used are inevitably somewhat noisy. 
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computations, since the multi-qubit case follows similar 
lines, and requires no new ideas. 

Consider a dangling node implementation of a single- 
qubit cluster-state computation. The literal circuit for 
such a computation is depicted in Figure 14*2*1 The first 
step of the circuit sets up the initial microcluster. The 
remaining steps of the circuit alternate between applica- 
tions of the operation ADD, which effectively adds two 
extra layers of qubits to the cluster, and the opera- 
tion MEASURE, which implements the desired single-qubit 
measurements and feedforward of measurement results. 
(Note that we have omitted ancillas, classical processing 
of data, and feedforward from the visual depiction, for 
simplicity.) 

The operation ADD may be broken down into the fol- 
lowing operations: 

• Preparation of microcluster states which we at- 
tempt to adjoin to the existing cluster. 

• All quantum gates applied in the process of at- 
tempting to adjoin the additional microclusters, in- 
cluding non-deterministic CPHASE gates. 

• Ancillas used in the non-deterministic CPHASE 
gates. 

• All the classical control and feedforward. 

• The computational basis measurements and local Z 
operations needed to remove undesired qubits from 
the cluster. 

In addition to these physical operations, to simplify 
the analysis it is convenient to append some fictitious 
perfect controlled-SWAP operations to ensure that the ef- 
fective output qubits always appear on the same output 
lines. That is, these controlled-SWAP operations effec- 
tively change the labels on the qubits to ensure that the 
same qubits always act as the output, but otherwise don't 
affect the output state. The operations are chosen so that 
the approximate circuit identity illustrated in Figure El 
holds. This identity illustrates the fact that the result 



FIG. 43: The operation ADD satisfies the approximate circuit 
identity illustrated here. This identity becomes more exact as 
k is increased, and the level of noise in the physical operations 
used to do the operation is decreased. 

The next step of the proof is to replace the measure- 
ments, classical control and feedforward performed in 
the ADD operation by equivalent ancilla preparations and 
unitary quantum operations, in a similar fashion to our 
threshold proof for deterministic CPHASE gates in Sec- 
tion ]W\ The result is an operation QADD, involving only 
quantum systems, perfect ancilla preparation, and noisy 
unitary operations. Note that, once again, we model 
noisy ancilla preparation by perfect ancilla preparation, 
followed by a noisy quantum memory step. The output 
of the noisy dangling node cluster-state computation is 
thus the same as the output of the imperfect circuit il- 
lustrated in Figure PHI 

As in the deterministic case, we may insert fictitious 
perfect extra gates between QADD operations without 









QADD 


MEAS. 






















QADD 


MEAS. 














i+>- 
[+)- 

|anc) -j- 




— r— 




-H 



FIG. 44: The output of this imperfect circuit is the same as 
the output of the noisy cluster-state computation. The state 
|anc) is an ancilla prepared (noisily) offline; it contains mi- 
croclusters, the ancillas used in the KLM CPHASE, and qubits 
used to simulate the classical control and feedforward in the 
dangling node implementation. Note that the ancilla line is 
marked by a slash to indicate that it involves many systems, 
not just a single qubit. 



25 



changing the overall output, as illustrated in Fig 05] The 
entire noisy cluster-state computation is thus equivalent 
to the repeating circuit illustrated in Figure where 
the gate G is as defined in Figure 1771 

In our analysis it is convenient to distinguish between 
the perfect gate G, and the imperfect gate G n , which also 
includes the effects of noise due to interactions with the 
environment, E. Note that due to the chaining property, 
Aq-.ae(G, G n ) < Nr), where the systems Q and A are as 
defined in Figure07J N is the total number of noisy oper- 
ations involved in the operation G n , and r\ is the maximal 
noise strength of any of those operations. Simple count- 
ing shows that N < c\k 2 for some positive constant c\ in 
the range 10°-10 2 , and so Aq : ae(G, G n ) < cik 2 rj. 

With these definitions we see that G is a unitary posts- 
elected gate acting on the registers Q and A, implement- 
ing the gate V defined in Figure 1381 on register Q, with 
probability p s = 1 — pl~ l . Let S be the subspace of 
states for Q of the form |^)|+) . . . |+), where is an 
arbitrary single-qubit state. Let T be the subspace of 
states for the system QA which are of the form |</>)|anc), 
where \<p) is any state in S, and |anc) is the initial state 
of the ancilla A. It is also useful to define T" to be the 
subspace of the combined system QAE containing states 
of the form \4>)se ® |anc), where \<P)se is an arbitrary 
state in S ® E. Applying Lemma ^ we conclude that 
there exists a unitary W acting on A such that 



\\G\ T -{V® W)\ T \\ = ^2(1- VPl)- 



(34) 



Recalling that Aq : ae(G, G n ) < cik 2 rj, we see that there 
exists unitary Ge acting on E such that 



\\G n -G®G E \\ < Cl k\ 
Restricting to the subspace X" we have 

\\G n \ T .-(G®G E )\T>\\ <crk\ 
From Eq. I|34[l we deduce that 



(35) 
(36) 



\\(G®G E )\ T . - {V®W®Ge)\t>\\ = yfltl 



(37) 
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FIG. 46: The output of this imperfect circuit is the same as 
the output of the noisy dangling node cluster-state computa- 
tion. 



Using the triangle inequality and Eqs. H3(jfl and (|37|l we 

obtain 



||G„| T ' - (V ® W ® G e )\t> II < c 1 k 2 V + ^2{l - ^j). 

(38) 

Applying the second unitary extension theorem, Theo- 
rem |21 we see that there exists unitary G n whose action 
on T' is identical to the action of G n , and such that 



\\G n -{V ®W ®G E )\\ < 2cifc"?7 + 2 x /2(f- VpI) 



(39) 



Summarizing, the output of the noisy cluster-state 
computation is the same as the output of the noisy cir- 
cuit in Figure The blocks of two imperfect CPHASE 
gates shown in this circuit represent the operation G n , 
which we have seen satisfies 



A Q:AE (V, G n ) < 2 Cl fc^ + 2 J 2(1 - VPT), (40) 



where p s — 1 — pjT is the probability of successfully 
adjoining a microcluster, and pf is the probability of a 
non-deterministic CPHASE failing. Examining Figure 1391 
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FIG. 45: We may insert some fictitious perfect gates between 
the operations in Figure B3I without changing the output. 



FIG. 47: The definition of the gate G used in FigureUjO This 
figure also defines subsystem labels Q and A. 
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— — f FigureEU However, the analysis for such operations goes 

through in exactly the same way as in the single-qubit 

case. 

— * — Summing up, we have shown that noise of strength 

up to 77 in a dangling node implementation of a cluster- 
state computation is equivalent to noise of strength at 
most cifc 2 7?+C2?7+2^/2(l — ^/pT) in the original quantum 

circuit, where c\ is a positive constant of order 10°-10 2 , 

C2 is a positive constant of order 10°, and p., = 1 — 2p k ~ l + 

FIG. 48: The definition of the gate V. Pf k ~ 2 '■ Thus, provided k and rj satisfy 



+>- 
+>- 
+>- 
+>- 



MEAS. 



cifc 2 ?7 + c 2 r] + 2 J 2(1 - VpT) < Vth 



(41) 



MEAS. 



anc 



l+>- 
1+)- 



anc 



1+) 

1+) 

I anc) — j- 



FIG. 49: The output of this noisy circuit is the same as the 
output of the noisy cluster-state computation. 



we see that a dangling node implementation using non- 
deterministic gates is equivalent to a noisy two-at-a-time 
deterministic implementation, where the additional lay- 
ers of the cluster are added with noise of total strength 



at most 2c\k 2 r\ + 2^/2(1 — ^JpT). Combining this with 
the results about the two-at-a-time implementation in 
Subsection IV CI we see that a noisy dangling node im- 
plementation of a single-qubit cluster-state computation 
is equivalent to a noisy single-qubit quantum circuit com- 
putation. Noise of strength rj in the dangling node im- 
plementation is mapped to equivalent noise of strength 
c\k 2 rj- 



-C2 77+2^/2(1 — s/Ps) in the quantum circuit model 
of computation, where the extra contribution cir\ is due 
to noise in the single-qubit measurements, and C2 is a 
positive constant of order 10°. 

Similar reasoning can be used to map noise of strength 
77 in a dangling node implementation of a multi-qubit 
cluster-state computation back to the original quantum 
circuit with noise of strength at most c\k 2 r\ + C2?7 + 
2^/2(1 — y/p2), where the constants ci and C2 may now 
be different, but are still of the same order, and p s — 
1 - 2p k f ~ 1 + p 2k ~ 2 . The only difference in the proof is 
that in addition to the G operations, multi-qubit com- 
putations also involve analogous operations based on the 
process of adjoining the more complex microcluster of 



it is possible to compute fault-tolerantly in the optical 
cluster-state proposal for computation. This may be 
rephrased as the condition: 



n < ti ocs 



77 th - 2^/2(1-^7) 



c\k 2 



01 



(42) 



We can always ensure that ry£ c t s h > by choosing k suffi- 
ciently large. Provided this condition is satisfied, ry^ c t s h is 
thus a threshold for optical cluster-state computation. 



VI. CONCLUSION 

In this paper we have proved two fault-tolerant thresh- 
old theorems for the cluster-state model of quantum com- 
putation. Our first threshold theorem applies to im- 
plementations in which deterministic (but noisy) entan- 
gling gates are available. Our second threshold theorem 
is specifically adapted to the case of optical quantum 
computation, where entangling gates are performed non- 
deterministically. In both cases our threshold theorems 
hold for quite pessimistic noise models, allowing non- 
Markovian noise applied by an intelligent adversary who 
can exploit constructive interference to enhance the ef- 
fects of the noise, and even cause errors in the classical 
computation and feedforward of measurement results. A 
drawback of our noise models is that they do not yet 
allow the possibility of leakage errors, like photon loss 
in optics. We expect to remove this drawback in future 
work, by combining the ideas in j2^ with well-established 
methods for dealing with photon loss, e.g. |42^ . 

Our focus has been on proving that a finite threshold 
exists for cluster-state computation, rather than on ob- 
taining a precise numerical evaluation of the threshold. 
This is consistent with our general philosophy of under- 
standing the threshold through a two-part process: first, 
rigorously proving the existence of a finite threshold for 
some large class of noise models; and second, through a 
combination of numerical and analytic work obtaining a 
realistic estimate of the threshold for some specific and 
physically-motivated noise model. In this paper we have 
obtained a rigorous proof that a finite threshold exists. 
Detailed numerical simulation and optimization of the 
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threshold value for realistic noise models is underway, 
and will be reported elsewhere [27|. 

Our investigations in this paper have been geared to- 
ward variants of the one-way quantum computing model 
introduced by Raussendorf and Briegel 1]. However, the 
techniques we have proposed seem quite generally appli- 
cable to the task of making measurement-based schemes 
for quantum computation fault-tolerant. It seems likely 
that schemes such as those proposed in 0, H, EJ El (see 
also references therein) can be made fault-tolerant us- 
ing similar ideas. Particularly appealing from a theo- 
retical point of view is the possibility of fault-tolerant 
computation using measurement alone, with no unitary 
gates whatsoever (excepting quantum memory). This 
might be done using, for example, a scheme such as 0, 
or one of the simpler variants that has since been pro- 
posed 00 0113 

The most important conclusion from our results is 
that noise need not be an obstacle to scalable quantum 
computation using cluster states. In particular, our re- 
sults provide encouraging evidence that practical propos- 
als for cluster-state quantum computation using neutral 
atoms and optics 01 are > m principle, fully scalable 
approaches to quantum computation. 
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APPENDIX A: THE UNITARY EXTENSION 
THEOREMS 

In this appendix we restate and prove the first and 
second unitary extension theorems, from Section III Al 

Theorem 5 (First unitary extension theorem) Let 

U, U and V be unitaries acting on a Hilbert space T . Sup- 
pose S is a subspace of T such that U and U have the 
same action on S, i.e., U\g = U\s- (Note that we do not 
assume that U and U leave the subspace S invariant, so 
U\s and U\s should be considered as maps from S into 
T .) Then there exists a unitary extension V of V\s to 
the entire space T such that 



It is worth noting that the proof below holds not just 
for the matrix norm, but for any norm such that ||^4-B| < 
||A|| and ||W|| < 1 for all unitaries W. 

Proof: Let P be the projector onto the subspace S, 
and let Q = I — P be the projector onto the orthocom- 
plement of S in T. Define 



V = VP + VU^UQ. 



(A2) 



We will show that V has the required properties. It is 
clear that V is an extension of V\s- To prove unitarity 
of V we first observe that 



rfUQ = QU^U. 



(A3) 



To prove this equation, observe that it is equivalent to 
WUP = PWtj, since Q = I - P. But WUP = PUW 
follows easily from the fact that U\s = U\s- The unitar- 
ity of V follows from Eqs. I|A2(I . (|A3|I . and some algebra: 



VPV^ + V{U^U)Q{U^U) V f 
VPV f + VQV ] 

I. 



(A4) 
(A5) 
(A6) 
(A7) 



To bound || — U\\ observe that 

V-U = VP + VU^UQ - UP -UQ. (A8) 
Observing that UP — UP and inserting UU' — I we get 
V-U = VP + VU^UQ -UP -UU f UQ (A9) 



{V -U){P + U^UQ). 



(A10) 



Recall that ||AB|| < \\A\\ \\B\\, and observe that P + 
WUQ is unitary. It follows that || P + U^UQ\\ = 1, and 
thus \\V- U\\ < ||V-I7||, as required. 
QED 

Theorem 6 (Second unitary extension theorem) 

Let U and V be unitary operations acting on a (finite- 
dimensional) inner product space T. Suppose S is a sub- 
space ofT. Then there exists a unitary operation V such 
that V\s = V\s and 



\U-V\\<2\\U\ s -V\s\\. 



(All) 



\\v-u\\<\\v-u\\. 



(Al) 



The proof of the second unitary extension theorem is 
somewhat more complex than the proof of the first. We 
begin the proof by introducing some notation related to 
the singular value decomposition of a matrix, and then 
we state and prove some simple lemmas about singular 
values and matrix norms. 

Recall that the singular value decomposition states 
that an arbitrary m x n matrix M can be written 
Al = Lm^mRm, where Lm is an m x m unitary ma- 
trix, Rm is an n x n unitary matrix, and Em is an 
to x n matrix, all of whose entries are zero, except the 
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diagonal entries (Ejf)jj = o~j(M), known as the singu- 
lar values, which are non-negative and arranged in de- 
creasing order, oi(M) > oi (M) > .... It is easy to see 
that the singular values are determined by the equation 
<jj (M) 2 = Xj(M^M) = XAMM*), the jth largest eigen- 
values of the matrices M T M and MM'. It is also useful 
to note that ||M|| = a\{M). When M is an m x m 
matrix (i.e., a square matrix), we will use the notation 
cr m i n (M) = er m (M) to denote the smallest singular value 
of M. 

Proposition 7 Suppose U is a unitary matrix that can 
be written 



U 



A C 
B D 



(A12) 



where A is an m x m matrix, B is n x m, C is m x n, 
and D is n x n. Then cr m i n (A) = a m i n (D). 

Proof: Inspection of the equations U'U = UU* = I 
implies that A* A + B*B = I m and BE* + DD* = I n , 
where we use Ik to denote the k x k identity matrix. 
Using these equations we have: 



= X m (A*A) 

= 1 - Xi(B*B) 

= l-Xi(BB*) 

= X n (DD*) 



(A13) 
(A14) 
(A15) 
(A16) 
(A17) 



from which the result follows. An alternate proof of this 
proposition follows immediately from the well-known CS 
decomposition of linear algebra (see, for example, Theo- 
rem VII. 1.6 on page 196 of |43j1. 
QED 

To state the next proposition we need to introduce 
some additional notation. We define the partial order 
X < Y for matrices X and Y if Y — X is a positive 
matrix. We define \X\ = VxTX. 

Proposition 8 Suppose M is a matrix such that \M\ < 
I. Then \\I - M\\ > 1 - cr min (Af). 

Proof: Choose a normalized vector \if>) so that 
Sm|V*) = a min(M)\'ip) . By the singular value decomposi- 
tion ||7 - M\\ = \\I - Lm^mRmW = \\W - where 
W = LljRij. Thus 

||/ -M|| > \\(W-^mM\\ 

> \\wm - \\z M \Tp)\\ 

= 1 - cr min (M), 

where we have used the triangle inequality. 
QED 



(A18) 
(A19) 
(A20) 



Proposition 9 Let 



M = 



A B 

£t c 



(A21) 



be a positive (and thus Hermitian) square matrix. We 
assume A is m x m, B is m x n, and C is n x n. Then 



W\\ < \\A\\ + \\C\\ 



(A22) 



Proof: This proposition can be viewed as a spe- 
cial case of Aronszajn's inequality, Theorem III. 2. 9 on 
page 64 of |4^. Alternately, since M is positive we can 
find a block matrix D = [D 1 D 2 ] such that M — D*D, 
and thus A = d\D\ and C = D ] 2 D 2 . We have 



M| 



= Xi(D*D) 

= \i(DDl) 

= X 1 (D l Dl-\ 

< Xi(DiD\) 

= \\A\\ + \\B\\ 



d 2 d\) 

-Xi{D 2 Dl) 



(A23) 
(A24) 
(A25) 
(A26) 
(A27) 



where the second-last line follows from the well-known 
eigenvalue inequality Xi(A + B) < X%(A) + X%(B), true 
for all Hermitian matrices A and B. 
QED 

Proof of the second unitary extension theorem: 

We prove the theorem in two parts. In the first part we 
prove the result for the case U = I. In the second part we 
show that the general result follows from the case when 
U = I. 

For the first part, we write V in block-diagonal form 

as 



V = 



A C 
B D 



(A28) 



where the first block represents a basis for 5 1 , and the 
second block represents a basis for the orthocomplement 
S±. Let the singular value decomposition of D be D — 
LdYidRd- Rotating by V D to change the basis of S± we 
see that V can be written in the new basis as 



V 



a' a 

B YidRdLd 



(A29) 



where A',B',C represent the action of V with respect 
to the new basis. We define the extension V by 



V 



V 



I 







o lJX 



A' C" 
B' S B 



(A30) 
(A31) 



where C" = C'L' D R' D . It is clear that V is unitary and 
V\s = Vs. All that remains to complete the first part of 
the proof is to bound \\V — I\\. We have 



\V-If = \\(V - I)(V* - I)\\ 

= \\2I-V-V*\\ 

21 -A 1 -A'* -C"-B* 
-B' - C"t 21 - 2S D 



(A32) 
(A33) 



(A34) 
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Applying Proposition [5] we obtain 

||y-/|| 2 < \\2I-A'-A'*\\ + \\2I-2E D \\. (A35) 

But \\2I-A'-A*\\ < \\I-A'\\ + \\I-A'm = 2||I-A'|| < 
2||7| s -y|s||. We also have ||27-2S D || = 2-2a min (D) = 
2-2a min (A') < 2\\I-A'\\ < 2\\I\ S - V\ s \\. The first part 
of the proof follows. 

For the second part of the proof we set U' = J, V = 
U'V. Then by the first part of the proof there exists 
unitary V' such that V'\ s = UW\ S , and || V"' - I\\ < 
2\\U^V\ S - Is\\- Set V = UV'. Then V is a unitary 
matrix such that V\s — V\s and 

\\V-U\\ = \\V'-I\\ (A36) 
< 2\\UW\ S - I\s\\ (A37) 
= 2\\V\s-U\ s \\. (A38) 

QED 

An alternate proof of the second unitary extension the- 
orem may be given, following similar lines, but based on 
the well-known CS decomposition from linear algebra. 
We have taken the approach presented here as it is only 
slightly more complex than the alternate proof, and relies 
on less background material. 



APPENDIX B: PROOF OF PROPOSITION [5] 

In this appendix we prove Proposition [5] which as- 
serts the truth of the circuit identity in Figure l2"3l 
This proposition establishes the correspondence between 
blocks QB a and gates HZ a in the quantum circuit. 

To prove Proposition[S]we make use of four simple cir- 
cuit identities which we assert here without proof. Each 
may be readily verified by direct computation. The first 
identity is a simplification of a common element in literal 
circuits for cluster-state computation. 



-©- 



-©■ 



-©- 



Identity 3 The following circuit identity holds: 



-e- 



-©■ 



= f 9 



Identity 4 The following circuit identity holds: 



-©- 



-©- 



= -ffi- 



With these identities we can now prove Proposition 
The proof is a short sequence of circuit identities which 
are shown in Figure 1501 



Identity 1 The following circuit simplification holds 
where the second qubit is in the state |+). 



~HZ n 



HZ n 



1+)- 



The remaining three identities concern the commuta- 
tivity properties of CNOT (controlled-NOT) and CPHASE 
gates when the control of one operator is the target of 
the operator immediately following or preceding it. 



Identity 2 The following circuit identity holds: 



9 T HZ ±c 



-©- 



By the definition of QB a in Figure 20. 



-f HZ± a ffi ©- 



-©- 



Using Identity 1. 



-ffi ©- 



-* — ©- 



Commute the cnot through the con- 
trolled HZ± a . 



~HZ a 



-e- 



Commuting the CPHASE through the 
HZ a and applying Identity 2 to the fi- 
nal tWO CNOTS. 



HZ a 



-e- 



Commuting the rightmost CPHASE 
through the two preceding CNOTs us- 
ing Identity 3. 



HZ a 



1+) 



-©- 



Commuting the rightmost CPHASE 
through the SWAP and cancelling with 
the leftmost cphase, followed by com- 
muting the remaining cphase leftward. 





HZ a 









1+) 















5 



Commuting either cnot through the 
swap and cancelling. This is our de- 
sired result. 



FIG. 50: The sequence of circuit identities used to prove Proposition |5] 
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